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ABSTRACT 

We investigate the metallicity evolution and metal content of the intergalactic medium 
(IGM) and galactic halo gas from z — 2 — > using 110-million particle cosmological 
hydrodynamic simulations. We focus on the detectability and physical properties of UV 
resonance metal-line absorbers observable with Hubble 's Cosmic Origins Spectrograph 
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(COS). We confirm that galactic superwind outflows are required to enrich the IGM 
to observed levels down to z = using three wind prescriptions contrasted to a 
no-wind simulation. Our favoured momentum-conserved wind prescription deposits 
CN| , metals closer to galaxies owing to its moderate energy input, while the more energetic 

■ constant wind model enriches the warm-hot IGM (WHIM) 6.4 times more. Despite 

these significant differences, all wind models produce metal-line statistics within a 
factor of two of existing observations. This is because O VI, C iv, Si iv, and Ne vin 
absorbers primarily arise from T < 10 5 K, photo-ionised gas that is enriched to similar 
levels in the three feedback schemes. O vi absorbers trace the diffuse phase with p/p < 
100, which is enriched to ~ 1/50 Z at z = 0, although the absorbers themselves 
usually exceed 0.3 Z and arise from inhomogeneously distributed, un-mixed winds. 
Turbulent broadening is required to match the observed equivalent width and column 
density statistics for O VI. C iv and Si iv absorbers trace primarily T ~ 10 4 K gas 
inside haloes (p/p > 100), although there appear to be too many Civ absorbers 
relative to observations. We predict COS will observe a population of Nevm photo- 
ionised absorbers tracing T < 10 5 K, p/p ~ 10 gas with equivalent widths of 10-20 mA. 
^ | Mgx and Sixn are rarely detected in COS S/N=30 simulated sight lines (dn/dz < 1), 

although simulated Sixn detections trace halo gas at T = 10 6 ~ 7 K. In general, the 
IGM is enriched in an outside-in manner, where wind-blown metals released at higher 
redshift reach lower overdensities, resulting in higher ionisation species tracing lower- 
density, older metals. At z = 0, the 90% of baryons outside of galaxies_are enriched 
to Z = 0.096 Z Q , but the 65% of unbound baryons in the IGM have Z = 0.018 Z 
and contain only 4% of all metals, a large decline from 20% at z — 2, because metals 
from early winds often re-accrete onto galaxies while later winds are less likely to 
escape their haloes. We emphasise that our results are sensitive to how metal mixing 
is treated in the simulations, and argue that the lack of mixing in our scheme may be 
the largest difference with other similar publications. 

Key words: galaxies: formation; intergalactic medium; cosmology: theory; methods: 
numerical 
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1 INTRODUCTION 

The intergalactic medium (IGM) contains the majority of 
the Universe's baryons. The best available IGM tracer is 
the intervening Lya absorption forest observed in quasar 



spectra |Rauchlll99sh . Contributing a lesser but still signif- 
icant amount of absorption are intervening me tal lines that 
trace the en r ichment history of the I GM ( e .g. iTvtleJ 1987 : 



ICowie et all 1 19951 : iBurles fc Tvtleil 1 19961 : ISongailal 



2001; 
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Schave et al. 120031; ISimcoe et all |2004| ; iDanforth fc Shulll 
20081 ; iTripp et all 120081 )" Since metal absorption is often 



found far from the sites of star formation, analytic models 
and simulations strongly suggest that metals must be 
transported out from galaxies i nto the IGM via large- 



19931; Nath & Trenthaml 19971; lAguirre et al.l 2001; 


Springel & Hernauist 


2003bl; lODDcnhcimer & Dave 


20061; Cen & Fane 


2006); iKobavashi et al. 


2007; 


Sommer-Larson & 


mbol 120081; Wiersma et al 


2009a; 


Tornatore et al.ll20ld 


Choi & Naeaminell201ll: Tescari et al. 


201l|). 



The observation and interpretation of IGM metal-line 
absorbers provides critical insights into the physics of galaxy 
formation. First, galactic outflows that enrich the IGM con- 
currently regulate the growth of galaxies by removing gas 
from star-forming regions. Second, the census of IGM metals 
provides a unique and independent constraint on the inte- 
grated amount of star formation in our Universe; indeed, the 
metal census within galaxies falls significantly short of the 
expectations from cosmic star formati on, which is known 
as the "missing metals problem" (e.g. iFukugita fc Peebles! 
l2004tlDave fc Oppenheimerjl2007r i. Finally, the amount and 
physi cal state of the metals in the IGM alters how gas cools 
fe.g. ISutherland fc Dopital 1 19931 ; iGnat fc Sternberd 120071 ; 
W lcrsma et al.l l2009b). thereby impacting the e fficiency by 
which galaxies receive gas from the IGM (e.g. Katz et al.l 
Il996t iKeres et all 2005; lOppenheimer et al. I l2010h . Hence 
studying IGM metal absorbers provides unique insights into 
the cycle of baryons between galaxies and their surrounding 
gas. 

Despite the importance of IGM metal absorbers, the 
chasm between observations and a theory of metal lines re- 
mains wide. Many basic questions concerning IGM chemi- 
cal enrichment remain poorly understood. Answering such 
questions is now becoming feasible owing to the recent in- 
stallation of the Cosmic Origins Spectrograph (COS) on the 
Hubble Space Telescope (HST), whose central mission is to 
improve our knowledge of the low-redshift IGM. Key ques- 
tions include: 

• What is the metallicity of the IGM? At redshifts of 
z = 2 — 3 measurements r ange between 10~ 3 ' 5 — 10~ 2 Zq 
dSongaila fc Cowiel Il99rj ; iDave et all 1 19981 ; ISchave et all 
120031 ; Isimcoe et alj|2004). while by z = a value of 0.1 Zq 
is typically quot ed (Cen et al.l 1200 ll ; iFang fc Brvanl l200ll ; 



IChen et aij|2003h . While this suggests an increasing metal- 
licity with time, the interpretation is complicated by an 
evolution in how metals trace the IGM baryon distribu- 
tion. Owing to cosmic expansion, a given column den- 
sity H i system traces a higher overdensity at lower red- 
shifts (|Dave et al.l [l999j); a similar trend, modulated by 
evolving photo-ionisation rates, also exists for metals. One 
consequence of these effects is that summing the column 
densities of metal lines to find an integrated ion density 
(fiion) does not necessarily translat e into a global metallic- 
ity fflmrtai). lOppenheimer fc Pavel |2006l ') argued that the 
non-trivial ionisation correction from C IV to total carbon 
abundance evolved with redshift, resulting in a relatively 
flat Q(j IV that masks a significant increase in the metallic- 
ity of the IGM. The se measurements have now been made 
down to z ~ (e.g. IDanforth fc Shulll |200S| ; ICooksev et al.l 
2010), and COS should provide significantly more insights 



into the relationship between metallicity and density over 
most of cosmic time. 

• What temperature is the enriched IGM gas? High- 
ionisation lines commonly observed in the IGM such as 
C IV and O VI can arise both from photo-ionisation and 
collisional ionisation. The low densities and hard quasar- 
dominated metagalactic flux can produce such high ion- 
isation states, as can hotter gas in the warm-hot inter- 
galactic medium (WHIM) of gas shock-heated during large- 
scale structure formation. There is vigorous debate as to 
whether observed high-ionisation lines trace collisionall> 
ionised hot metals at T > 10 5 K (e.g. 
Danforth fc Shulll |200S|; ICen fc Chisaril 



Cen fc Fang 2006 



20111 ISmith et all 



201 it iTepper-Garcia et all l201lj ~ or low-density photo 
ionised gas at T < 10 s K (e.g. iThom fc Chenl l2008bl ; 
ITripp et all 120081 lOppenheimer fc Pavel 120091 . hereafter 
OP09). The increased sensitivity of COS enables multiple 
ions to be detected in many systems, allowing independent 
constraints on the nature of high-ionisation systems. 
• How homogeneously are metals distributed in the IGM? 
If metals are well-mixed throughout the IGM, metal lines 
will cool more baryons, possibly leading to increased gas 
accretio n onto galaxies and h ence an increase in star for- 
mation (W lcrsma et al.ll2009al). If metals are poorly mixed 
as ob servations suggest ( Simcoe et al.l 120061 : ISchave et alj 
12007ft . then a census of O VI cannot st raightforwardly be 
converted to a census of baryons (e.g. ITripp et al l l200d ; 
IDanforth fc Shulll [2005 ) since Ovi only traces the baryon s 
in over-enriched regions (OD09; iTepper -Garcia et al . 2011). 
The analysis of observed metal-line systems often yields 
inconclusive answers, with comple x absorbers suggesting 



multi-phase gas components (e.g. P rochaska ct al. 2004 



2006; Co oksev et ail 



[ Savage et all 120051; Ibanforth et al . . 
boost bowk et al.ll2009l ; lNaravanan et al.ll2009l ). The higher 
signal-to-noise provided by COS will be critical for disen- 
tangling these effects. 

A number of reasons exist for our less complete knowl- 
edge of metal absorbers compared to that of the Lya for- 
est. Observationally, metal lines are weaker and, therefore, 
harder to detect. Moreover, many key ionisation states are 
buried within the Lya forest or occur at X-ray wavelengths 
where current sensitivities are limit ed (e g lChen et al.l l2003l; 
iNicastro et al. [ 120051 ; ICen fc Fangj 120061 : IFang et all l201fj| ). 
Theoretically, IGM enrichment depends sensitively on the 
transport mechanism out of galaxies, which is poorly con- 
strained, while t he Lya forest is l ess affected (Kollm eier et 
al. 2003; OD09; IDave et al.ll2010tlTescari et al.ll20ilh . Fur- 
thermore, converting metal ions into metallicities typically 
requires a detailed knowledge of the photo-ionisation rates 
at energies poorly constrained by the data or, if collision- 
ally ionised, is high ly sensitive to gas temperatures (e.g. 



Heckman et al. 2002) and hence detailed cooling rates (e.g. 
Gnat fc Sternberg! 120071 ). However, there are many more 
metal transitions observable compared to hydrogen transi- 
tions, providing an opportunity to synthesise a wider suite of 
observations into a single coherent framework to determine 
the physical state of IGM metals. Hydrodynamic simulations 
that evolve the metallicity and ionisation state of the IGM 
self-consistently along with the enriching galaxy population 
offer a way to develop such a framework within a modern 
hierarchical context. 

In this paper, we employ state-of-the-art cosmological 
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hydrodynamic simulations to study the z — — 2 IGM as 
traced by high ionisation metal lines. Our strategy involves 
examining what can be learned by observations targeting 
quasar absorption line spectra randomly sampling the IGM, 
which represents a volume-weighted measure. We focus on 
high ionisation lines including O VI, C IV, Si IV, Ne VIII, Mg x, 
and Si XII. Low ionisation lines (e.g. C II, Si II, Mg n) are 
not addressed as they trace higher density gas where self- 
shielding renders the assumption of a uniform ionisation 
background insufficient. We briefly comment on N v a nd O V 
predictions. This paper directly follows Dave et al. |2010l . 
Paper I in this series, hereafter D10), which discussed the 
Lya forest and the physical conditions of baryons in the 
IGM. The next paper in this series will study correlations 
between galaxies and absorbers (Kollmeier et al. in prep., 
Paper III). 

While we will show predictions for observable quanti- 
ties with COS, we do not intend these to be used for direct 
comparisons with data. We are more interested in showing 
how plausible variations in the physics manifest themselves 
in forefront observations, particularly in terms of answering 
some of the questions outlined above. For instance, we will 
discuss how integrated column densities trace the underlying 
IGM metallicity, how various ions observable with COS trace 
the temperature and density of the enriched IGM baryons, 
and how absorber alignment statistics can constrain the ho- 
mogeneity of metals in the IGM. We will do all this for 
several prescriptions of galactic outflows, to better under- 
stand how COS observations can constrain this important 
physical process of galaxy formation. 

A description of our cosmological hydrodynamic simula- 
tions and our artificial spectra follows in §2. §3 concentrates 
on the evolving physical conditions of metals, subdivided 
by phase between z — 2 — > 0. We simulate COS metal-line 
absorber observations in §4 and relate them to their phys- 
ical conditions in §5, demonstrating that only a fraction of 
IGM metals are probed via UV resonance absorption lines. 
We spend §6 exploring physically plausible model variations 
both to demonstrate the inherent uncertainties in interpret- 
ing metal-line observations and to distill what can be learned 
using relatively straight-forward statistics in a sample size 
similar to the one that COS will observe. §7 discusses the 
successes and failures of our simulations and compares them 
to simulations by other groups. We conclude in §8. For read- 
ers who want a briefer overview of the results, we suggest 
first reading the conclusions section and following the sec- 
tions listed after each of the main findings for more detail. 



2 METHODS AND SIMULATIONS 

We utilise our modified version (|Oppenheimer fc Davel2008T ) 
of the Gadget-2 N-body + Smo othed Particle Hydrody- 
namic (SPH) code l|S P ringelll2005l l to evolve a series of cos- 
mological simulations to z = 0. The details of the simula- 
tions can be fo und in §2 of DIP; they hav e also been used 
and described in lOppenheimer et al.l(|2010h . B riefly, our cos- 
mology agrees with the WMAP-7 constraints jjarosik et al.l 
2011) and uses the following values: SIm = 0.28, JIa = 0.72, 
h = Ho /(100 kms _1 Mpc _1 ) = 0.7, a primordial power spec- 
trum index n = 0.96, an amplitude of the mass fluctuations 
scaled to a 8 = 0.82, and Q b = 0.046. 



2.1 Wind Implementations 

We briefly recap the implementations of galactic winds 
in these simulations; f or a m ore extensive description see 
lOppenheimer fc Pavel j2008). The basic technical imple- 
mentat ion is the one described in ISpringel fc Hernquistl 
l|2003al ) and incorporated into Gadget-2. In a given 
timestep, a particle in a star-forming region has p e j — r\ x 
SFR/m p of being ejected, where r\ is the mass-loading fac- 
tor, SFR the star-formation rate, and m p the particle mass. 
The ejection velocity, « w i n d, is oriented according to v x V$, 
where v is the particle velocity and V<3> is the gradient of the 
gravitational potential, resulting in winds propagating per- 
pendicular to galactic disks. To allow the particle to escape 
the ISM of its host, hydrodynamic forces are turned off until 
one of two criteria are met: either 1.95 x 10 10 kms _1 / u wiiid 
years have passed or (more often) the particle flows out to a 
density of less than 10% of the star formation critical den- 
sity. 

In our constant wind (cw) model, we employ an ejec- 
tion velocity of u W md = 680 kms -1 and a mass loading 
factor of r\ = 2 for all galaxies. The slow wind model 
has i> w i n d = 340 km s" 1 and r? = 2. Both of the se wind 
models are similar to ISpringel fc Hernquistl (|2003al ) . except 
for their assumed velocity; they use « w ind = 484 kms -1 . 
The momentum-conserved wind mo del, inspired by the 
iMurrav. Quatert. fc Thompson! (|2005l ) models of radiation 
pressure-driven winds, has « W md = 3(7 \fJL — 1 where /l is 
the luminosity in units of the Eddington luminosity required 
to expel gas from a galaxy potential, and r\ — ao/a where 
(To = 150 kms -1 and a is the galaxy's internal velocity dis- 
persion calculated by a group finder run during the evolution 
of the simulation (see lOppenheimer fc Davej|2008T ). 

The metallicity of an ejected wind particle is the metal- 
licity of the SPH particle at the time of ejection. SPH par- 
ticles are continually self-enriched while they are eligible 
for star formation. In this scheme, both ejected winds and 
newly formed stars have the metallicity of the SPH par- 
ticle at the time of the ejection or formation event; wind 
particles are essentially ejected ISM SPH particles. Note 
that only particles above the critical star-formation density 
threshold are eligible to form stars or to become launched 
as a wind. Cooling, which is never suppressed, invariably 
allows the launched winds to cool and maintain a temper- 
ature of ~ 10 4 K while decoupled. Once the hydrodynamic 
forces are restored, which always occurs deep inside haloes 
given the two criteria listed above, the SPH particle inter- 
acts normally. In practice, re-coupled wind particles show 
a bimodal behaviour: (i) they shock heat to temperatures 
above the peak cooling efficiencies for metal-enriched gas 
(T ~ 10 5 — 10 6 K), or (ii) they shock heat to a lower temper- 
ature where metal-line cooling is efficient and then rapidly 
cool to T ~ 10 4 K. The most important determining factor 
appears to be w w i n d, which we will show when we consider 
IGM enrichment patterns among the different wind models 
in S 



2.2 Simulations 

Our naming convention for our simulations continues 
the precedent we established in our recent papers: 
r[boxsizc]n [number of particles per side] [wind model], where 
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the initial letter "r" indicates the particular choice of cos- 
mology above. The main simulations that we explore are 
evolved in a cubic periodic volume that is 48/i _1 Mpc on a 
side: r48n384nw (nw: no winds), r48n384cw (cw: constant 
winds), r48n384sw (sw: slow winds), and r48n384vzw (vzw: 
moment um-conserved winds ) . The se simulation s were also 
used in Oppenheimer et al.l l|2010T i. D10, and iDave et al.l 
The main simulations are listed in Table [1] The 
gas and dark matter particle masses are 3.56 x 1O 7 M and 
1.81 x 10 s Mrs,, respectively. 

Four other simulations listed in Table [1] are evolved in 
a 16/i _1 Mpc box with 2 x 128 3 particles, the same reso- 
lution as the main simulations. We use these to test the 
effects of metal-line coolin g with photo-ionisation on m etal- 
line observa t ions, because Tepper-Garcia et al.l (|201ll ) and 
I Smith et alj ()201 lh demonstrate that this could significantly 
alter observations of Ovi versus using collisionally ionised 
metal-line cooling rates. Despite their smaller volume, we 
will demonstrate in £16.21 that these simulations' results re- 
main converged compared to the 48 ft. _1 Mpc simulations, 
likely because absorbers are mainly associated with sub-M* 
galaxies that are sufficiently sampled in both size volumes. 
We add the suffix PI to those simu lation names evolved usi ng 
the photo-ionised cooling rates of IWiersma et all (|2009bh . 

2.3 Spectra 

Spectra are generated as described in §2.3 of D10. Briefly, 
70 continuous lines of sight are created from z = 2 — > for 
each model using our spectr al generation code specexbin 
|Oppenheimer fe Pavel I2006T I. We generate mock spectra 
using only the strongest transition for an ion (e.g. the 
1031.93A line for Ovi, the 1548.19A line for Civ) and con- 
volve the spectra with the COS line spread function (LSF). 
We add Gaussian random noise with a signal to noise ratio 
S/N=30 per 6 kms -1 pixel to simulate the highest qual- 
ity data obtainable using COS. We employ the continuum 
provided directly by the simulation, and do not attempt to 
mimic the observational procedure of fitting a continuum, 
since at low redshifts the sparseness of absorption features 
usually results in a robust and reliable continuum fit. 

We fit absorbers in the spectra using Au- 
to VP l|Dave et al l Il997n , obtaining a column density 
TV, line width b, and redshift z for each absorber. We 
consider each metal transition separately and do not com- 
bine them into a single spectrum with all the transitions. 
This assumes that, in observations, each transition can be 
reliably identified and fit without contamination from other 
lines, which is typically true at low redshifts. 



2.4 Ionisation Background 

The metagalactic flux amplitude is adjusted to match the 
observed evolution of the Lya forest mean flux decrement. 
We multiply the Lya optical depth by the same values of 
/ T that we calculated in D10 and presented in Table [1] 
of that paper. In p ractice, we renormalise our CLOUDY 
IIFerland et al. I ll998l . version 08.00) ionisation tables so that 
the magnitude of the ionisation background is multiplied by 
1// T , which effectively shifts the ionisation fractions of all 
ions by 1// T in density under th e assumption of an opti- 
cally thin gas. The vzw normalised Ifiaardt fc Madaul |200ll . 
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Figure 1. The ioni sation background spectr a explored in this pa- 
per. HM2001 is the lHaardt fc MadaJ j200ll) spectrum multiplied 
by 1.5 to match the observed Lya flux decre ment (DIP). Three 
other backgrounds generated using CUBA llHaardt fc Madaul 
l200lt) with varying slopes for quasar spectra above the Lyman 
limit are also displayed (HM2005). Energies are listed in units 
of Rydbergs as well as electron volts (eV, top x-axis). Below are 
photo-ionisation cross-sections for various ions. Coloured regions 
show the magnitude of the cross-section from the ionisation state 
below to the specified ionisation state, while black regions show 
the cross-section from the specified state to the state above. For 
these panels, levels of 10 — 18 cm 2 are indicated by large ticks on 
the j/-axis and 10 — "cm 2 are indicated by small ticks, with the 
i/-axis starting at zero in all cases. 



HM2001) z = 0.2 background (i.e. J„ multiplied by 1.5) is 
shown in Figure[T]in black. Three other backgrounds are also 
displayed, which will be described in H6. II when we explore 
the effects of variations in the ionisation background. 

We also display the photo-i onisation cross-sections of 
various ions used in CLOUDY l|Verner fc Yakovlevl Il995l : 
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Table 1. Simulations 



Name- L b v wind d rf E wind /E f SN Z Cooling^ U. t /£l b (z = 0)* 



r48n384nw 


48 


384 






0.0 


CIE 


0.9 


0.213 


r48n384cw 


48 


384 


680 


2 


0.96 


CIE 


0.9 


0.045 


r48n384sw 


48 


384 


340 


2 


0.24 


CIE 


0.8 


0.124 


r48n384vzw 


48 


384 


oc a 


OC (T -1 


0.45-0.56 


CIE 


0.67 


0.097 


rl6nl28cw 


16 


128 


680 


2 


0.96 


CIE 


0.9 


0.040 


rl6nl28cw-PI 


16 


128 


680 


2 


0.96 


PI 


0.9 


0.037 


rl6nl28vzw 


16 


128 


OC cr 


OC CT _1 


0.41-0.54 


CIE 


0.67 


0.096 


rl6nl28vzw-PI 


16 


128 


oc a 


OC <7 _1 


0.41-0.56 


PI 


0.67 


0.088 



"Suffix: nw- no- wind, cw- constant wind, sw- slow wind, vzw- momentum-conserved winds. 

b Box length of cubic volume, in comoving h — 1 Mpc. 

c Numbcr of dark matter and SPH particles per box side. 

d Wind velocity in kms -1 . 

e Mass loading, where M w ; n[ j = rjM*. 

^Wind energy divided by total SNe energy assuming Esn = 9-7 X 10 48 erg/ Mg. Range of values at z =0,1,2,3 for vzw. 
ff Metal-line cooling. CIE- Collisional Ionising Equilibrium cooling rates, PI- Photo-ionisation-dependent cooling rates. 
h Factor Lya optical depth is scaled in determining ionisation background normalisation. 
' Fraction of baryons in stars at z = 0. 



IVerner et all 1 1996ft on a linear scale in the bottom panels 
of Figure [T] The coloured regions show the ionisation cross- 
section of the next lowest ionisation state to the specified 
ionisation state, while the black shading indicates the cross- 
section from that state to the next highest state. For exam- 
ple, the magenta shading indicates the cross-section for Si ill 
to Si IV, and the black shading indicates the cross-section for 
Si IV to Si v. With the exception of Si IV, all cases shown are 
lithium-like ions and display a trend of increasing ionisa- 
tion energy and declining cross-section with higher atomic 
number. Metal ion cross-sections are complex, which is ex- 
emplified by the Si in and Si IV cross-sections peaking above 
10 Ryd even though Sim has the lowest ionisation energy 
(2.5 Ryd) of all the species we consider. The multiple offset 
wedges, two in the case of lithium-like ions, correspond to 
different subshells in the ground state. 

2.5 Chemodynamical Model 

We follow the production of four heavy elements (C, O, Si, & 
Fe) from three sources: Type II Supe rnovae (SNe), Type la 
SNe, and AGB stars as explained in lOppenheimer fe Pavel 
(2008). The details of the enrichment scheme are in 62.1 of 
D10. Cosmic oxygen and silicon production is dominated 
by Type II SNe, while carbon production is augmented 
signifi cantly by AGB stars but s till primarily arises from 
SNe (|Qppenheimer fe Pavel feoOS). We do not follow neon 
and magnesium directly when generating Nevin and Mgx 
statistics but instead use oxygen as their proxy, which is 
a reasonable approximation given that all three elements 
arise primarily from Type II SNe. The re sulting abundances 
are_[Ne/O]=0.16 and [Mg/Ol =-0.21 for IChieffi fe Limongl 
(2004) SNe yields when using lAsplund et al. I |2005l ) solar 
abundances, which we scale to throughout this paper when- 
ever quoting metallicities. We assume the same Type II SNe 
origin for nitrogen by using [N/O]=-0.40, but note that this 
does not include secondary nitrogen production in AGB 
stars, which is likely significant. The solar metal mass frac- 
tion by weight, Zq, is 0.0122 and the mass fractions of 



carbon, nitrogen, oxygen, neon, magnesium, and silicon are 
0.00218, 0.00062, 0.00541, 0.00103, 0.00061, and 0.00067 re- 
spectively. We also do not consider depletion of these ele- 
ments onto dust grains, althou gh recent results suggest that 
there may be dust in the IGM ^Menard et al.ll20ld ; TZu et all 
l201lT ) that contains as much as 50% of the mass of intergalac- 
tic metals. 



3 METAL PHYSICAL CONDITIONS 

3.1 Density- Temperature Phase Space 

The density-temperature phase space diagrams of metals re- 
siding in the IGM are displayed in the left panels of Figure [2] 
for our preferred vzw model at z = 2, 1, and 0. The diagrams 
are coloured by metallicity while the brightness indicates the 
fraction of baryons. On the right side are depicted the z = 
phase space diagrams for the other three wind models. His- 
tograms along the x- and y-axes depict the relative fraction 
of metal mass at a given density and temperature, respec- 
tively. 

The metals populate the cosmic phase diagram differ- 
ently than the baryons (see e.g. P10) with the density and 
temperature dependence of metallicity linked to the type of 
outflows that put them there, which we will discuss in H3.2I 
Here we focus on generic trends seen among all the wind 
models, contrasting them with the no-wind case. 

Following P10, we divide the gas component into four 
phases, demarcated by the solid lines in the figures: 

• Piffuse (T <T th ,S < 5 th ), 

• WHIM (T>T tb ,8 <S th ), 

• Hot halo (T > Tth, 5 > 5 t h), 

• Condensed (T < T th , S > <5 t h), 

where T t h = 10 5 K represents the division between the 
"cool" and "hot" phases, and 5th represents the division 
between bound (halo) gas and unbound (IGM) gas. These 
definitions are designed to provide a physically-motivated 
demarcation between cosmic baryon phases. Tth = 10 5 K 
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is near the peak in the helium and metal cooling curve 
that causes a minimum in the temperature distribution of 
baryons, while the density threshold separates truly inter- 
galactic baryons, i.e. those not associated with gravitation- 
ally collapsed structures, from baryons that lie predomi- 
nantly within the viria l radii of galaxy halos. Following 
IKitavama fc Sutol l| 19961 ). we set 

5 th = 6tt 2 (1 + 0.4093(l//n - l) ' 9052 ) - 1, (1) 
where 



J On(l + z) 3 + (1 - fi m - Q A )(1 + Z) 2 + OV ' W 

5th evolves from a value of « 60 at high redshift to « 120 
by z = 0. As in DIP , we a dopt Equation [1] to be l/3rd 
the IKitavama fc Sutol l| 19961 ) mean enclosed density, which 
corresponds to the density at 7-200 for objects just undergoing 
collapse. The condensed phase does not include star-forming 
ISM gas, which has nu ^ 0.13 cm -3 and lies to the right of 
the high-density edge of each diagram. The dashed vertical 
line indicates the mean overdensity of the Universe (S = 
p/p — 1 = 0). The percentages listed provide a fractional 
accounting of the total metals in all gas and star particles 
associated with each phase; we also show the percentage of 
metals in galaxies in the lower right corner, and additionally 
divide this phase into stars and and star-forming (ISM) gas. 

Examining the evolutionary trends in the vzw model, 
we see that the fraction of metals in the diffuse IGM drops 
rapidly with time, from 18% to 3% from z = 2 — ¥ 0. In 
contrast, the fraction in the WHIM phase is quite small 
(~ 1 — 2%) and evolves downwards more slowly. Recall from 
D10 that the mass fraction of the WHIM in this model in- 
creases from 10 — > 24% during this time, while the fraction 
in the diffuse gas drops from 73 — ¥ 41%. Hence overall, the 
fraction of metals in the IGM (i.e. diffuse+WHIM) drops, 
a trend that is accentuated in the diffuse phase by a drop 
in the baryonic mass fraction in this phase, while for the 
WHIM the trend is countered by an increasing mass frac- 
tion. The decrease in the IGM metal fraction is balanced by 
an increase in the metal fraction in galaxies, from half at 
z ~ 2 to almost 80% at z = 0. Hence, from z = 2 — > 0, frac- 
tionally the metals move from low-density regions to higher 
density regions. We emphasise, however, that galaxies con- 
tinue to enrich the IGM up to the present day and that the 
metallicity of the IGM does not actually decline; only the 
fraction of metals in the IGM declines. We return to this 
point in H3.3I 

By z — 0, we see that all the wind models have most 
of their metals within galaxies. In the no-wind case, only 
3% of all metals are outside of galaxies, while in the cw 
case it is 37%. The momentum-conserved and slow wind 
models are intermediate, with 10-20% of their metals out- 
side of galaxies. Still, these z — numbers are not large; 
in all the wind models more metals are outside of galax- 
ies at high- z. Thi s means that the solution to the "missing 
meta ls" problem (|Bouche et al.ll2007l ; lDave fc Oppenheimerl 
120071 ) at early epochs depends more on the accounting of 
metals outside of galaxies than today. 

Examining the temperature histograms along the y-axis 
(including only metals outside galaxies), we see that the gas- 
phase cosmic metals in all the models at all epochs have a 
maximum at T ~ 10 4 K. This is an artificial temperature 



floor, since we do not include radiative cooling below 10 4 
K; however, this also is nearly the equilibrium temperature 
where cooling equals photo-heating. We will show in H6.2I 
that when we exa mine models incl u ding m etal-line cooling 
below 10 4 K from IWiersma et alj (|2009bl ). metals cluster 
around ~ 10 4 K, although with more variable temperatures 
extending above and below owing to different equilibrium 
temperatures dependent on density and metallicity. 

The central result is that most cosmic metals are not 
in hot gas, either in the WHIM or hot halo gas. At z = 0, 
the IGM metallicity distribution for all wind models has a 
bimodal temperature distribution, with cool metals mostly 
at T ~ 10 4 K and hot metals mainly residing at intraclus- 
ter medium (ICM) temperatures (T > 10 6 K). Few met- 
als reside in between 10 5 K and 10 6 K (except in the cw 
case) because metal-line cooling is very efficient there, caus- 
ing metal enriched gas to rapidly transition through these 
temperatures (OD09). In £17.21 we state that the bimodal- 
ity is stronger in our simulations compared to simulations 
that explicitly smooth metals, because metals confined to 
individual SPH particles cool more rapidly over this inter- 
mediate temperature range. 

The density histograms along the bottom of Figure[2]for 
the vzw model (left panels) show a peak in metal fraction at 
nu ~ 10~ 3 — 10~ 4 cm" 3 at all redshifts. This corresponds to 
an increasingly higher o verdensity towards lower reds hifts, 
reflecting the result from lDave fc Oppenheimer] (2007) that 
the peak overdensity of metals increases with time. This 
is qualitatively true in the other wind models as well (not 
shown), although the actual distribution of metals in over- 
density depends strongly on the particular wind model as 
we discuss in £13.21 

Figure [3] shows the mass fraction (left panels) and 
mass-weighted mean metallicity (right panels) of six bary- 
onic phases (stars-top, ISM-2nd, diffuse-3rd, WHIM-4th, 
hot halo-5th, & condensed-bottom panels) for our four wind 
models, from 2 = 4 — > 0. At early epochs, star-forming gas 
holds the majority of metals. As time passes, metals become 
progressively locked up into stars, so that by 2 = stars con- 
tain the majority of cosmic metals. Outflows suppress the 
metal content in stars and the ISM somewhat but do not 
alter the overall trends. The mean stellar and ISM metal- 
licities increase slowly from 2 = 4 — > 0, but even at early 
epochs it is already close to solar, and by today it exceeds 
solar. This is consistent with the early, rapid enrichment 
in galaxies typical in these models (e.g. iDave et al.1 120061 ; 
iFinlator etldl feoiO) and is relevant for the IGM because 
it shows that the material expelled from galaxies has near- 
solar metallicities at all epochs. The metallicity of the IGM 
depends then on how far the metals are driven, and how 
they mix with the unenriched ambient gas. 

The metal fraction of the IGM phases (diffuse+WHIM) 
generally drops from 2 ~ 3 —> in the wind models. This 
reflects the increasing difficulty that metals have in escap- 
ing from the denser regions where they are formed to less 
dense regions, as depicted in Figure [2] While the fraction 
of metals in the IGM phases is generally not large, it is 
still sufficient to reproduce the observations of metal-line ab- 
sorbers across cosmic time fe.g. [Oppenheimer fc Dave1l2006l ; 
lOppenheimer. Dave, fc Finlatorll2009l )~ 

For the WHIM and hot halo phases, the 5 t h threshold 
that divides these two phases, first introduced in D10, makes 
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Figure 2. Dcnsity-tempcraturc phase space diagrams showing the distribution of metals divided into four phases by the solid lines: 
horizontal divides "hot" and "cold" phases, and vertical divides halo and IGM phases. Colour indicates metallicity, and brightness 
corresponds to the baryon fraction. Percentages indicate the fraction of global metals in each phase (bottom right- condensed, bottom 
left- diffuse, top left- WHIM, top right- hot haloes, plus the percentage in galaxies subdivided into stars and ISM is also listed). Grey 
histograms show the gaseous metal fraction as a function of density (x-axis) and temperature (y-axis) including only metals outside 
galaxies. The vertical dashed line corresponds to the mean overdensity of the Universe. The left panels show the evolution of the 
momentum-conserved wind (vzw) model at z = 2, 1, and the right panels show the three other wind models at z = (nw- no winds, 
sw- slow winds, & cw- constant winds). 



a stark difference in the WHIM metal content compared to 
previous WHIM definitions without such a density thresh- 
old. The WHIM, which contains 24% of z = vzw baryons, 
holds only 1.2% of the z — metals and remains steady 
at [Z/H]~ -1.8 since z = 4. The WHIM forms through 
shock-heating owi ng to the gravitationa l collapse of large- 
scale structures (ICen fc Qstrikerl 1 19991 ; iDave et all 1 19991 . 
l200ll ; ICen fc Qstrikerl 120061 ) while metal-enriching galactic 



outflow s contribute fractional l y less to the WHIM at lat e 
times (|Oppenheimer fc Pavel 120081 ; I Cen fc Chisaril l201ll ). 
The metal content of the hot halo phase steadily increases, 
and by z — it holds more metals than the combined IGM 
phases (cf. 5.2% vs. 3.9% for vzw). The y-axis histograms for 
the vzw model clearly show the growing importance of hot 
halo metals from z — 2 — > 0. The z = panel indicates that 
most of these metals are at T > 10 6 K and hence are in the 
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Figure 3. The evolution of metal mass fractions (left) and mass- 
weighted mean metallicity (right) of the four wind models from 
z = 4 — + subdivided by phase (lst-stars, 2nd-ISM, 3rd-diffuse, 
4th- WHIM, 5th- hot halo, k. 6th-condensed). Four wind mod- 
els are shown: solid blue- no-wind, dashed green- momentum- 
conserved winds, dotted red- constant winds, and long-dashed 
orange- slow winds. 

intragroup or intracluster medium (ICM). In all the wind 
models, the hot halo metallicity reaches 0.1 Zq at z = 0, 
but it also shows a density gradient such that the portion 
that would emit most strongly in X-rays has [Z/H]~ —0.5 
(see Figure [2]). This i s comparable to the measured metal- 
licity of the ICM (e.g. iFukazawa et al.|["l998l ; iPeterson et al.l 
120031 ) and reconfirms that these simulatio ns can generally r e- 
produce the observed ICM metallicities (|Dave et al.| [2008). 

More broadly, the fraction of metals in the ICM in- 
creases at very early epochs as winds enrich diffuse cos- 
mic gas (|Oppenheimer. Dave, fc Finlatorll2009l '). but then at 
later epochs the metals do not escape their haloes so easily 



or else re-accrete back into haloes. Hence, the metal fraction 
in the diffuse phase peaks at z ~ 3 — 4 in these wind mod- 
els (Figure [3| . Meanwhile, the metallicity in all the phases 
increases, but slowly; the most rapid increase is in the con- 
densed phase gas, since this halo gas continues to be enriched 
even to late times by outflows. These trends, while gener- 
ally true, are in certain aspects quite sensitive to the wind 
model, as we will discuss next. 

3.2 Dependence on the Galactic Wind Model 

The value of studying enrichment patterns lies in their high 
sensitivity to galactic winds, upon w hich the process of 
galaxy formation heavily depends (e.g. ISchave et al.l l20ld : 
lOppenheimer et alj|201fj) . We now contrast our vzw simu- 
lation to the three other wind models. The no-wind model 
produces nearly twice as many stars and metals as any wind 
model, and is incapable of enriching the diffuse IGM appre- 
ciably at any redshift, leaving it mostly pristine even by 
z = (Figure [2] upper right). The necessity to enrich the 
diffuse IGM to levels of at least [Z/H] = -3 at z > 2 (e.g. 
ISchave et al.l 120031 ; ISimcoe et al.l 12004) is the main empir- 
ical requirement for large-scale galactic outflows. The dif- 
fuse phase metallicity of the nw model remains below the 
y-axis in Figure [3] only reaching [Z/H] = —6.1 at z = 4 
and —3.9 by z = 0. This model does, however, achieve a hot 
halo metallicity comparable to our other feedback models, 
indicating that the enrichment of group and cluster haloes 
at least in part owes to a mechanism other than galactic 
outflows, such as dynamical stripping. I Dave et al.l (2008) 
showed that their no-wind model reached observationally 
reasonable [Fe/H] abundances via delayed Type la SNe 
from intragroup/intracluster stars, but [O/H] abundances 
were too low compared to observations and the baryonic 
fraction in stars was too high. Overall, the no-wind model 
stron gly fails to match basic observations of both galaxies 
(e.g. iDave et al.ll2011ah and the IGM. 

The constant wind (cw) model has a high wind speed 
of fwind = 680 kms -1 emanating from every galaxy, which 
pushes shock-heated metals into voids. This leads to a 
WHIM metal mass fraction more than 10 times higher than 
any other model (Figure middle panel). These metals re- 
main hot, because adiabatic cooling from Hubble expansion 
is not sufficient to return this gas to the photo-ionised IGM 
equa tion of state (EoS; iHui fc G ncdin 199 J ISchave et al.l 
2000) within a Hubble time. The excess heating is also re- 
flected in the phase space diagram for this wind model (Fig- 
ure [51 lower right), which has a high- metallicity peak in un- 
derdense regions at T ~ 10 6 K. The fast winds also in- 
crease the WHIM baryon percentage from 23 to 33% (com- 
pared to the nw model), whereas the vzw and sw mod- 
els increase the WHIM bary onic fraction by one tenth a s 
much, a 1% increase (D10). As lOppenheimer fc Davi (|2006h 
show, a similar model with somewhat lower wind speeds 
(■y W md = 484 kms -1 ) produces Civ absorbers at z ~ 3 that 
are are too wide compared to the observations, indicating 
that there is too much IGM heating at this early epoch. 

The slow wind (sw) model with v w ind = 340 kms -1 does 
not overly heat the early IGM, and in some sense represents 
a hybrid of the vzw model (similar average wind velocities) 
and the cw model (same mass loading). Like cw, this model 
creates a concentration of shock-heated metals (T ~ 10 4 ' 5 
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Figure 4. Mass-weighted mean metallicity as a function of over- 
density from z = 3 — > for the four models: momentum-conserved 
winds (vz w) , no winds (nw) , constant winds (cw) , and slow winds 
(sw). 

K) in voids, but it is at lower temperatures and is much 
less prominent owing to the sw model's lower wind speeds. 
In this model low-mass galaxies eject winds at high veloci- 
ties relative to their escape velocity and high-mass galaxies, 
typically in dense environments, eject winds at less than the 
escape velocity. Hence, the sw model cannot populate the 
ICM with metals and, therefore, the metal phase space dia- 
gram in this regime mimics that of the no-wind model. 

In the vzw model, winds from small galaxies have lower 
velocities. Hence they deposit metals at higher densities 
where metal-line cooling is efficient and can return these 



metals to T ~ 10 4 K. In hot gas, the vzw model produces 
a strong metallicity-density gradient by z — owing to its 
high outflow velocities from massive galaxies. Winds with 
> 1000 kms" 1 speeds emanating from these galaxies rarely 
escape the halo owing to the intervening halo gas, but they 
do shock heat to temperatures high enough that the cooling 
times are long, and hence the metals do not re-accrete onto 
the galaxies. 

Figure |3] shows the mean metallicity-density relation in 
our four wind models from z — 3 — > 0. The no-wind model 
fails to provide widespread IGM enrichment, but even dy- 
namical processes can deposit a smattering of metals at the 
cosmic mean density by z ~ 0. All the wind models show 
increased metallicities at all the plotted overdensities rela- 
tive to the no-wind case. The constant wind model results 
in a curious upturn in the mean metallicity at void densi- 
ties (p < p), which is a generic feature of strong feedback 
models (e.g. ICen fc Chisaril I20TH . their Figure 24). Unfor- 
tunately this is very difficult to probe observationally using 
UV resonance lines owing to the low densities, as we will 
discuss in SJS] The slow wind model has a less pronounced 
upturn, and it shows a less widespread metal distribution 
particularly at early epochs owing to its lower wind speed. 
The vzw model shows an interesting trend with overden- 
sity, and the metallicity in filaments (1 < p/p < 100) is 
nearly time-invariant. Continual enrichment from galaxies 
pumps metals into these regions, but this is approximately 
balanced by gravitational effects that draw the metals back 
closer to galaxies and by the accretion of more pristine 
gas. At higher overdensities, the mean metallicity increases 
with time, while at sub-mean densities, metals s lowly diffuse 
deeper into voids l|Oppenheimer fc DavejlioOfj ). The differ- 
ences between the models offer an opportunity to probe out- 
flows if the metallicity-density relation and its evolution can 
be reliably measured, which is a central goal of IGM enrich- 
ment observations. 

Overall, the evolution of metal content by phase shows 
much more dependence on feedback than our corresponding 
study of the abundance of baryons by phase (D10). While all 
the wind models significantly suppress the global star for- 
mation rate and hence metal production, the distribution 
of metals in the IGM shows clear differences. The velocity 
of the outflow and its motion through the ambient gas de- 
termines to what temperature the gas is shock heated and 
to what overdensity it reaches. Gas in sufficiently dense re- 
gions cools quickly, particularly since it is enriched. Hence, 
the outflow parameters critically determine the enrichment 
history of baryons outside of galaxies, providing an oppor- 
tunity to constrain such parameters based on observations 
of IGM metal lines. 

3.3 Outside-In IGM Enrichment 

A generic feature of our IGM enrichment models is that 
they have a very early epoch where metals are distributed 
outwards into the IGM, but over most of cosmic time the 
newly produced metals became confined to ever decreasing 
regions around the galaxies that produced them. We call this 
"outside-in" IGM enrichment, and it has consequences not 
just for how metals enrich the IGM, but for how wind m a- 
terial re-accretes onto galaxies l|Oppenheimer et al.ll2010h . 
We have already seen trends exemplifying outside-in 
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enrichment above. The metal mass fraction in the IGM 
decreases steadily since z ~ 3, and the peak overdensit 
of metals increases with time. lOppenheimer fc Pavel (1200, 
showed that the median distance reached by vzw outflows 
(with a wide dispersion) is ~ fOO physical kpc, implying 
widespread enrichment (in a comoving sense) early on and 
more confined enrichment at later epochs. The evolution 
of the vzw metallicity-density relationship in Figure [4] (top 
panel) shows that metals are in place at IGM overdensities 
at 2 = 3 and that the metallicity is mostly invariant and 
actually declines slightly at p/p ~ 2 — 10 3 from z = 2 — > 0. 
Newly produced metals are mostly confined to their own 
haloes and fall back onto their parent galaxies primarily 
from these overdensities; however winds continue to enrich 
IGM densities until z = 0, especially the winds from lower 
mass galaxies, often even reaching void densities. Addition- 
ally, some enriched WHIM material also falls into galaxy 
haloes and eventually onto the central galaxy. 

Another view of outside-in enrichment is provided by 
the age of z = metals as a function of overdensity, shown 
in Figure top panel. We define the metal age as the time 
since a wind particle has most recently been launched, which 
is a good indicator of the age for metals in the IGM as these 
baryons are almost entirely enriched by winds. The median 
age of wind-blown metals is shown by the solid lines, and 
the 16 and 84% cuts are indicated by the dashed lines for 
the vzw model. 

Figure [5] shows directly that lower overdensities are en- 
riched earlier, particularly in the vzw model. For this model, 
the typical age of metals in voids is ~ 10 Gyr, reaching a 
minimum of ~ 1 Gyr at overdensities of a few thousand. 
The constant-?) models show less dependence on age with 
overdensity in the IGM, though at higher overdensities the 
metal ages are again significantly younger. The evolution of 
the metallicity-density relation of the cw model (Figure [4} 
shows a significantly increasing metallicity within voids with 
time. Most metals were already in place at higher redshift 
in gas that continues to expand into voids by z = 0, but 
galaxies continue to enrich void overdensities below z = 1. 
The same trends are true for the sw model, but they are 
shifted to higher overdensities and younger ages owing to its 
weaker winds. While the median age trends with overden- 
sity are statistically significant and indicative of the history 
of metal enrichment, the larger scatter in all models at most 
densities prohibits a direct mapping of an overdensity to an 
age for a given wind mode l. 

|w lersma et al.l (|2010h find a similar age-density anti- 
correlation using their definition of mean enrichment red- 
shift. Our definition of age is a lower limit compared to 
theirs, because we measure the most recent wind launch 
as opposed to the redshift-weighted enrichment age. Both 
methods will find similar ages for wind particles that rapidly 
enrich and launch once, as is typical for metals that escape 
galaxies at high-z, but we will find significantly younger ages 
for winds that recycle back onto galaxies and have multiple 
wind launches and enrichment epochs. 

To relate the IGM metal age to observables, we plot in 
the lower panel of Figure [5] histograms of the z = vzw cos- 
mic ion densities (in 0.1 dex bins of overdensity) for the ions 
that we follow and for H I. In general, higher ionisation metal 
species trace older IGM metals occupying lower overdensi- 
ties in the vzw model. For example, observed Nevill and 
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Figure 5. The top panel shows median ages of wind-blown metals 
at z = as a function of density (solid lines) with the 16th and 
84th percentage cuts (dashed lines, vzw only) for the momentum- 
conserved winds (vzw, green), constant winds (cw, red), and slow 
winds (sw, orange). Lower overdensities are enriched earlier in 
all wind models in what we term "outside-in" enrichment. The 
bottom panel shows cosmic ion densities of vzw metal species 
as a function of density, as well as, Hi (purple lines) and total 
metal density (grey lines). The two panels can be compared to 
show that in general higher ionisation species generally trace lower 
overdensities, which have older metals. 



Ovi absorbers are likely to be tracing fairly old (> 5 Gyr) 
metals, while C IV and Si IV more often trace younger metals, 
although with a large age scatter. Every metal species shows 
a peak corresponding to the overdensity where the species is 
primarily photo-ionised. A second peak at higher overdensi- 
ties for Si xil and Mgx corresponds to collisionally ionised 
metals inside haloes. Ne viii also shows an extended shoulder 
toward halo densities corresponding to collisionally ionised 
metals. This figure foreshadows our analysis in ^5] showing 
that species including Civ, Ovi, and Nevill are primarily 
photo-ionised in our simulations. Interestingly, these three 
species exceed the cosmic density of H I at a range of IGM 
densities. 

The idea that the IGM is enriched in an outside-in 
mann er has a long his to ry ( Tcgmark, Silk , fc Evrardl 
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Scannapieco et al. 2002T) . In general, these models postu- 
lated that the IGM underwent widespread early enrichment 
typically during the reionisation epoch from some putative 
powerful source such as Population III stars, when halo 
potential wells were still small. Our model, while similar 
in character, is fundamentally different in that it does 
not invoke very early enrichment from mysterious sources. 
Instead, the outside-in pattern arises from the dynamics 
of galactic outflows from the epoch of early galaxies 
(z > 6) until today. The source of IGM enrichment is the 
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observed population of high-z star-for ming galaxie s that 
our simulations n aturally reproduce (|Dave et al.l 120061; 
iBouche et aD 120071 ) . Energetically moderate outflows from 
these galaxies enrich a small, but rapidly growing fraction 
of the IGM, in agreeme nt with observed constrai nts from 
high-z metal-lines (e.g. iRvan- Weber et all I2009T ). These 
early outflows end up enriching a greater comoving volume 
than lower redshift outflo ws, which travel similar ph ysical 
distances in our model (lOppenheimer &: Pavel I2008T ). In 
contrast, lOppenheimer et al.l (|2007T l showed that enriching 
the IGM entirely prior to z ~ 6 would violate observations 
of the C IV enrichment at those epochs. 



4 METAL-LINE OBSERVABLES 

We present metal-line predictions for COS observations in 
this section along with comparisons to other pre-COS data. 
We consider statistics assuming a S/N=30 per 6 kms -1 
pixel in our simulated spectra, convolved with the COS line 
spread function (LSF) assuming the FUV G130M LSF at 
1450A, just as we did in D10. We model the LSF as a cen- 
tral peak with a la Gaussian width of ~ 17 kms" 1 and 
substantial non-Gaussian wings, keeping this fixed for all 
wavelengths to avoid introducing artificial evolution from a 
varying LSF. 

We obtain line statistics by fitti ng Voigt profiles t o indi- 
vidual components using AutoVP l|Pave et al.lll997t ). Each 
component has a fitted rest equivalent width (W), a column 
density (N), and linewidths (6-parameters). An absorption 
system is defined as all components that can be grouped to- 
gether if they lie within 100 kms -1 of another component. 
Systems are more easily comparable with observations, be- 
cause data quality does not always allow the easy and con- 
sistent identification of components. For example, it can be 
easier to divide a system into more components if the S/N 
and resolution are higher. We sum up system column den- 
sities to calculate cosmic ion densities (fi). We discuss the 
equivalent width distribution of systems versus components 
for O VI, where observations for both have been quantified. 
We present results in this section for our four main simu- 
lations, plus a vzw model with turbulent broadening added 
after the completion of the simulation (OP09), motivated 
by the observed O VI linewidths. 
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4.1 Equivalent Width Distributions 

Cumulative equivalent width distributions (EWPs), pre- 
sented in Figure [(J show the frequency of systems greater 
than or equal to the W indicated on the x-axis. They offer 
a simple and direct statistic to test theoretical models. We 
show EWPs of four metal lines: Ovi 1032A, Civ 1548A, 
Si IV 1394A (all between z = 0-0.5), and Nevm 770A (at 
z — 0.5 — 1.0). We show predictions for our four wind mod- 
els and a fifth model (vzw-turb in cyan) of turbulent 6- 
parameters (6 tur b) added after the completion of th e sim- 
ulation using specexbin l|Oppenheimer fc~P avc 2009). 

The EWPs with no winds (blue solid lines) exhibit far 
fewer absorbers than models with winds, and significantly 
underestimate the observations of Ovi, Civ, and Siiv. In 
general, the cw model (red dotted lines) produce fewer ab- 
sorbers than the other wind models because this wind model 



Figure 6. Cumulative equivalent width distributions of Ovi 
(top), Civ (second), Siiv (third), and Nevm (bottom) for 
our four wind models (solid blue- no-wind, dashed green- 
momentum-conserved winds, dotted red- constant winds, and 
long-dashed orange- slow winds, plus the vzw model with tur- 
bulent broadening (dot-dashed cyan). The line types correspond 
to t he same models throughout the paper. Observations for 
Ovi iDanforth fc Shulj|200g| ; llripp et aljESooih and Civ & Siiv 
llCooksev et alj|20ld . 1201 ll) are also plotted. For the first three 
panels, we count all lines in our continuous synthetic spectra out 
to z = 0.5 and divide by Az = 0.5, while the fourth panel uses 
the same procedure for z = 0.5 — 1. Error bars correspond to 1-a 
Poisson errors in the models and observations. 
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(i) produces less metals overall owing to its greater suppres- 
sion of star formation and (ii) injects more of its metals into 
the low-density WHIM where they are unobservable using 
these transitions. Compared to vzw, the sw model creates 
fewer Civ and Siiv absorbers, which we will show in [J5] 
trace the condensed phase within haloes that are enriched 
to higher levels by vzw winds. At the same time there exist 
more weak Ovi and Nevill lines in the sw model, because 
these arise from the diffuse phase that is enriched to higher 
levels via sw winds (cf. vzw and swz = metallicity-density 
relations in Figure [If. 

We discuss turbulence as a model variation in H6.3I but 
we present it here because this is part of our favoured model 
for Ovi observations. We also show the effect on Civ, Siiv, 
and Nevill, although we argue later that this model may 
not be applicable to absorbers tracing higher densities. As 
in OD09, adding foturb < 60 kms" 1 (see Equations [4] & [5] 
in H6.3P increases the number of absorbers above 100 mA , 
allowing us to fit the observed O VI EWDs, and it also blends 
some weaker components into stronger ones, thus reducing 
their numbers (cf. vzw and vzw-turb model in the top panel 
of Figure [6]l . Civ and Siiv also show a dramatic increase in 
the number of ^ 100 mA components if we apply this same 
density-dependent model. 

The simulations presented here produce statistically 
identical O VI results to OD09 for the vzw model. How- 
ever, the vzw model appears to over-produce the observed 
C IV, while the cw m odel provides a better fit to this data 
|Cooksev et all bOld fl. We predict that COS will make 
many Nevill detections, especially between 10-20 mA. In 
general, the lithium-like ions are harder to detect for higher 
atomic numbers corresponding to higher ionisation species, 
because both their oscillator strengths and wavelengths de- 
cline, leading to smaller equivalent widths for a given column 
density. 



4.2 Column Density Distributions 

Column densities provide a more direct accounting of the 
amount of metals present than equivalent widths, but they 
are more difficult to obtain from observations as COS often 
does not resolve individual metal-line absorbers. Neverthe- 
less, Auto VP fits column densities N and linewidths b for 
every component. Given that we have applied COS instru- 
mental characteristics (albeit with a perhaps optimistic S/N 
ratio), these predictions should be representative of the level 
of information obtainable from (the very best) COS spectra. 
We emphasise that these are distributions of inferred column 
densities obtained by Voigt-profile fitting components and 
summing them together into systems, and not necessarily 
distributions of true column densities of physically identi- 
fied, unblended components. 

We display the differential column density distributions 
(CDDs) of systems in Figure[7] As in D10, we plot f(N ioa ) = 
d 2 n/dNi OI1 dz multiplied by JVi on to improve the readability 
and to provide a dimensionless number that indicates the 



1 Kathy Cooksey kindly provided redshift frequency 
distribution fits for C IV and Si IV, which can be 
found at |http://w ww. ucolick.org/ ~xavicr/HSTCIV/ and 
http://www.ucolick.org/~xavier/HSTSiIV/ 



relative number of lines per unit redshift for a given Ni on . 
Again our four basic wind models are plotted, along with 
turbulent broadening added to the vzw model. We show 
eight species over Az = 0.5 for a redshift range that falls 
within the COS far-UV (FUV) channel: Ovi, Civ, Siiv, & 
Nv all at z = - 0.5, Nevill at z = 0.5 - 1.0, and Ov, 
Mgx, & Sixn at z = 1.0 - 1.5. 

The qualitative differences between the wind models 
and the data for the column density distributions are sim- 
ilar to those obtained for the equivalent widths as we dis- 
cussed in the previous section. Without winds, very few IGM 
metal lines are present, and the predictions fall well short 
of availabl e observations. The O V I CDD is consistent with 
the data l|Danforth fc Shulll l200ct FI for the vzw and vzw- 
turb models, but the same models produce too many C IV 
absorbers, especially above 10 14 cm -2 , compared to the ob- 
servations l|Cooksev et al.l l2010h . This disagreement is sta- 
tistically significant, and further data from COS should pro- 
vide much better constraints on weak C IV absorbers. Si IV 
predicts more strong lines and less weak lines than those 
observed (jCooksev et aljfeoi ll ). though the discrepancy is 
not as large as with C IV, and the = 0—1.2 redshift range 
of this data set is larger and centred at much higher red- 
shift than in our simulations ({z) ~ 0.9 vs. 0.25). Despite 
the different redshifts for Siiv, our simulated Siiv CDDs re- 
main statistically identical with (z) = 0.75 and 1.25 over the 
column density range of detected absorbers (not shown). 

We show Nv CDDs but caution that we do not fol- 
low the unique chemodynamic origin of nitrogen from sec- 
ondary enrichment sources (AGB stars). Hence, the Type 
II SNe origin assumed here for nitrogen likely represents 
a lower lim it, as n i troge n is copiously produced in AGB 



e.jecta (e.g. Marigol 200 ll) . Our models underestimate the 
iDanforth fc Shulll (l2008r ) statistics for higher Nv columns, 
which appears to support this interpretation. 

For Nevill (upper right panel of Figure 0, pushing 
below 10 13 ' 8 cm~ 2 is rewarding because the CDD rises 
sharply. At column densities > 10 14 cm -2 , the line counts 
are predicted to be well under one per unit redshift, but 
at ~ 10 13 cm -2 there should be 5 — 10 per unit redshift. 
Hence high S/N is essential to detect significant numbers 
of these absorbers. Note that our estimates are probably 
optimistic because we assume detection based on only the 
stronger component, but in actuality detection of the weaker 
780A component is usually required for a robust line iden- 
tification. 

We predict Ov should be generously abundant in the 
spectra of z > 1 quasars probed by COS over a wide range 
of column densities. We show the 630 A Ov statistics, be- 
cause it is likely representative of a large number of lower 
ionisation species with far-UV transitions that should con- 
tribute significantly to intervening absorption in higher red- 
shift quasars probed by COS. Ov should be among the 
strongest owing to a large oscillator strength and oxygen 
being the most abundant metal. 

Mg x and Si XII absorbers are very rare in our sight lines 
and definitely require winds to be present at all. We find 
only 2 Mg x and 6 Si XII absorbers over Az = 35 with W 



2 Charles Danforth kindly provided slightly updated versions of 
their data for O VI, C IV, Si IV, and N v. 
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Figure 7. Column density distributions for our mock S/N=30 COS observations for six species: Ovi 1032A, Civ 1548A, Siiv 1394A, 
Nv 1239A, Nevni 770A, Ov 630A, Mgx 609A, and SiXII 499A. The four wind models and the turbulently broadened vzw model are 
generated from a path length of Az = 35 o ver a redshift range observable by the COS NUV channel. Data from lDanforth fc Shull fcOOSl ) 
(triangles) and lCooksev et alj {2010, 2011) (dashed lines, thin lines encompass the 1-cr error range) are also shown. Our simulations do 
not include AGB production of Nv, so these predictions should be taken as lower bounds. Error bars are 1-cr Poisson errors. 



10 mA in the vzw model. We argue in i]5.1l that we may 
not be resolving the conductive interfaces from which very 
high-ionisation lines such as Nevni and Sixn may arise, and 
hence that these absorbers may be more common than our 
current models predict. 

Looking at the vzw-turb model, in theory the total col- 
umn density should be conserved (unlike W) when adding 
turbulence in the way that we do using specexbin. But as 



one can see, in practice the column densities are increased 
when turbulence is added. This occurs because for strong 
saturated lines, which preferentially receive the largest tur- 
bulent contributions, AutoVP tends to fit a higher column 
density given only a single line with no doublet informa- 
tion. Also, weaker lines may become blended or harder to 
detect. Hence, there is a noticeable increase in the predicted 
line frequencies especially for Civ at 10 140 cm -2 and Si IV 
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at 10 13 ' 8 cm -2 , which correspond to the column densities 
where these transitions fall off the linear curve of growth 
and begin to become saturated without turbulent broaden- 
ing. COS spectra will be able to use the two doublet tran- 
sitions to perform a curve of growth analysis and obtain 
more accurate iVs for saturated components. Including such 
information in AutoVP may mitigate the increase seen at 
high-iV; we leave such an analysis for a more careful com- 
parison to real data. 



4.3 Cosmic Ion Densities 

Integrating over a CDD or summing the column densities of 
lines provides a global view of the amount of cosmic metal 
absorption. To obtain the cosmic density of an ion, we sum 
up all the absorber column densities between 10 12-15 cm" 2 
using 



ion, z — 



H m mctal £iV(ion, z) 
cpcrit AX(z) 



(3) 



where Ho is the Hubble constant, c is the speed of 
light, pcrit is the critical density of the Universe, and 
^species is the atomic weight of the given species. AX(z) 
is the pathlength in a ACDM Universe where X(z) — 

^/f2ivi(l + z) s + J7a. We evaluate Equation [3] from the 
Voigt-profile fitted column densities, just as is done for ob- 
servations. We ch oose the colum n density range both for 
historical reasons (|Songailall200lh . and because it quantifies 
lower-column metal lines that are expected to arise in the 
IGM as opposed to those within the galactic ISM. 

We bin the f2i on s in bins of Az — 0.5 and plot the 
points with error bars in Figure [8] for Ovi, Civ, Siiv, and 
Nevm. For simplicity, we use the same COS FUV instru- 
ment parameters for the entire 2 = — 2 range even though 
most redshifts fall outside of the COS wavelength range. We 
show the same four wind models as before, and also the 
vzw model with turbulent broadening. We also plot f2 me tai 
computed by summing up the total mass density of all the 
SPH particles outside of galaxies (tih < 0.1 cm" 3 ) for the 
appropriate metal species (continuous lines). 

f2i on s usually do not change by much if we adjust the 
column density bounds from 10 12 to 10 13 cm -2 or from 
10 15 cm" 2 to any higher value, because the vast majority 
of an ion's cosmic density resides in 10 13-15 cm -3 compo- 
nents in S/N=30 spectra. The most dramatic exception is 
ClV at z > 1, which can increase by a factor of more than 
two when including lines above Nq iy > 10 15 cm" 2 . 

Models with outflows produce roughly constant f^Qvi 
and Q,(j IV from z = 2 — > 0, and agree with observations 
where available. While O VI agrees with the data in all three 
observables thus far explored, Civ only agrees for ^q iv 
because too many weak and too few strong absorbers ex- 
ist compared to the current data l|Cooksev et a.1.1 [20101 ) . In 
general, the non-evolution in these species reflects the non- 
evolution in the total mass densities of metals outside of 
galaxies (f2 me tais, lines in Figure [HJ, but the ion densities 
are shifted lower by a factor of ~ 30 — 100. The majority of 
the Universe's metals are produced after 2 = 2, but this is 
offset by more metals remaining nearer to galaxies and being 
reincorporated into stars, i.e. "outside-in" enrichment. Note 
that the ratio Sl mota i / fiion cannot be treated as a global ioni- 
sation correction for an ion (e.g. C/C IV or O/O Vl), because 
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Figure 8. Cosmic ion density evolution of Ovi (top), Civ (sec- 
ond), Si IV (third), and Nevm (bottom) for the same mod- 
els described in Figures \E\ and [7] Large symbols show simu- 
lated O; on s calculated from S/N=30 COS spectra, with hori- 
zontal error bars indicating Az = 0.5 bins and vertical error 
bars indicating the variance among Az = 0.5 sight lines. Con- 
tinuous lines show the evolution of f2 mcta i calculated by sum- 
ming up all the mass outside galaxies of oxygen (top panel), 
carbo n (second) , silicon (third ), and neon (bott om). Data are 
from iDanforth fc Shunl <200Sl) (solid triangle s ). ICooksev et al.l 
d2O10L l201ll) f solid squa res), fp'Odorico et al.| ll2010h (solid cir- 
cles). [ Pettini et al.l ||2003| ) (open triangle), and iRvan- Webe r et al.l 
(2009) (open square) and are all scaled to a Qjm = 0.28, Q/± = 
0.72, Hq = 70 kms" 1 cosmology. Redshifts are listed on the top 
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we will show in £|5.2I that f2j on calculated using Equation [3] 
does not recover the total cosmic density of that ion. 

In contrast to the other ions shown, ^Naym shows a 
steady decline from z = 1.5 — > for the vzw model. The 
photo-ionised Ne viii component is within the detection lim- 
its of COS at higher redshift, and becomes harder to de- 
tect at low redshifts at wavelengths only accessible using 
the Space Telescope Imaging Spectrograph (STIS) and the 
Far UV Spectroscopic Explorer (FUSE). Observing Ne vm at 
wavelengths between I100-1700A with COS could be ideal 
to trace the peak in ^N e ynp 

The turbulent vzw model predicts higher Qs than the 
vzw model even though the total column density theoret- 
ically should be preserved. As explained previously, the 
broader profiles with turbulent broadening are fitted more 
accurately by AutoVP than unbroadened thin saturated 
profiles. This is especially true if the cosmic density of an 
ion is dominated by a few strong absorbers. 

In summary, outflows are generically required to en- 
rich the IGM as observed, and there is a mild sensitivity 
in observations of W, N, and f2j on to the particular form 
of galactic outflows employed. While momentum-conserved 
winds plus turbulence best matches the O VI data (as it was 
in part designed to do), it overproduces Civ and Si IV ab- 
sorption. Nevm shows a particularly steep column density 
distribution owing to a preponderance of weak photo-ionised 
lines, motivating deep observations to detect this interest- 
ing IGM tracer. The cosmic metal density is generally con- 
stant from z ~ 2 — > 0, in contrast to it s mark ed rise from 
2 ~ 6 -> 2 ijQppenheimer fc Pavel 120061 . 120081 ). and this is 
broadly reflected in the constancy of Qi on for most ions over 
this redshift range. 



4.4 Component Linewidths 

Component linewidths (6-parameters) are also fit with Au- 
to VP, but are almost always dominated by the COS LSF 
(the instrumental b ~ 24 kms -1 ); therefore, linewidths are 
often much wider than they are in STIS data. We show 
the linewidth histograms for the simulations with winds in 
Figure [9] along with the vzw-turb model, which uses the tur- 
bulent broadening model motivated by OD09 to match the 
observed &-parameters of STIS Ovi absorbers, but now ap- 
plied to all ions. The histograms apply to all Ovi, Civ, and 
Si IV absorbers with W ^ 30 mA between z — — 0.5, and 
Nevm absorbers with W ^ 10 mA at z — 0.5 — 1.0. 

Without turbulent broadening, the linewidths of the 
three wind models are nearly identical. The COS LSF makes 
it difficult to determine temperatures from thermal broad- 
ening in our case since we show in the next section that the 
majority of our absorbers are near T ~ 10 4 K, correspond- 
ing to linewidths much smaller than the LSF. We do not 
show the nw model owing to its low number statistics. 

Adding turbulence using the heuristic prescription of 
OD09 dramatically alters the 6-parameter histograms for 
not just O VI, to which this prescription was calibrated, 
but even more so for C IV and Si IV. AutoVP fits a max- 
imum b — 50 kms" 1 , which means that the majority of 
Civ and Siiv absorbers in the 45-50 kms - bin may even 
be wider. Unlike thermal broadening, turbulent broaden- 
ing is independent of the atomic species, and the broader 
linewidths for C IV and Si IV are a result of the OD09 turbu- 
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Figure 9. Simulated COS linewidth histograms for Ovi (top), 
Civ (second), and Si IV (third) between z = — 0.5 for all com- 
ponents with W ^ 30 mA and Nevm (bottom) for components 
between z = 0.5 — 1.0 and W ^ 10 mA. The three models with 
winds and the turbulently broadened vzw model are generated 
from a path length of Az = 35 over a redshift range observable 
by the COS NUV channel. Note that these linewidths are not cor- 
rected for the instrumental broadening by the COS LSF, which 
has b Rd 24 kms -1 . 
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Figure 10. Cumulative equivalent width distributions of Ovi 
components and systems quantified in o ur simulation s for th e vzw 
and vzw-turb models. Observations of iTripp et al.1 J200ST) have 
both system (filled square) and component ( open square) distri- 
butions. Data from [Banforth fc Shulll <j200ST l (triangles) groups 
components based on the associated H I identification. 



lence model, which adds wider bturb to higher density gas. 
We argue in £|6.3l that such wide absorbers likely represent a 
break-down of the turbulence model, because while there are 



line p rofiles are apparent in the data (e.g. ICooksev et al.l 



few published observed linewidths of these s pecies, thinner 
e p r 

l20ld ). Ne vm absorbers are not as sensitive to the turbu- 
lence model, because they arise from lower densities where 
feturb is small. 



4.5 Components &: Systems 

The choice to group components into systems or treat them 
individually is an important consideration when comparing 
simulations to observations. K. Cooksey (private communi- 
cation) advises us that it is more proper to compar e grouped 
systems to the data of ICooksev et al.l l|20ld . l201ll ). because 
they group com ponents into systems if th ey lie within a few 
tens of kms - MDanforth fc Shull ( 20081 ) groups metal-line 
components together into systems if their associated Lycv 
lines lie within 50 k ms -1 of each other (C. Danforth, pri- 
vate communication). ITripp et al.l ()2008h publish both com- 
ponent and system EWDs, which we show in Figure [lOl 

The vzw model indicates how important it is to com- 
pare the appropriate component and system EWDs. The 
vzw system EWD would agree with the component EWD of 
ITripp et alJ l|2008l ). but the appropriate comparison, the vzw 
component EWD, under-predicts the strongest observed 
O VI components. The shortfall of strong W components was 
part of the motivation for OD09 to attempt the vzw-turb 
model (the other motivation being mismatched linewidths). 
The vzw-turb makes less of a difference between the two 
types of EWDs, because turbulent broadening increases the 
component Ws so that the strongest component of a sys- 
tem becomes more dominant relative to weaker components; 
thus the strongest component more likely sets the system W 
resulting in less differences between the two methods. The 
trend o f this model's comp onent and system EWDs matches 
that of ITripp et al.1 (|2008h . although we appear to consis- 
tently over-estimate component and system frequencies at 
«S 100 mA. 



The comparison to lDanforth fc Shulll l|2008r i is more dif- 
ficult because their search relies on the identification of the 
associated Hi absorption, but th eir measurements appear 
consistent with lTripp et al. I (|200Sl ) systems. We also tested 
the choice of grouping components within 250 kms -1 in- 
stead of 100 km s _1 , but this makes almost no difference and 
gives statistically identical results. If components lie close to 
each other, they are almost always within 100 kms -1 of an- 
other component. 



5 METAL-LINE ABSORBER PHYSICAL 
CONDITIONS 

We now link the observed metal-line absorber character- 
istics with their physical conditions just as we did for 
Lya absorbers in §5 of D10. We concentrate on ion mass 
density-weighted physical parameters, mainly hydrogen den- 
sity and temperature, using the same method described 
in OD09. This weighting procedure sums the ion mass 
density, and hence absorption, of every SPH particle con- 
tributing at the line centre of an identified absorber, ac- 
counting for peculiar velocities and temperature broadening. 
This procedure is similar t o the optical depth-weighting of 
iTepper-Garcia et al.l (|201ll ) with the only difference being 
that they sum physical quantities from each pixel across a 
line profile whereas we us e the value at the line centroid. 
ITepper-Garcia et al l (|201ll ) find that the difference between 
these two methods is small, so the methods should be com- 
parable. We consider all absorbers where the strongest line 
transition has W ^ 10 mA in our 70 S/N=30 COS sight 
lines. 



5.1 Absorbers in Phase Space 

Figure[TT]displays total metals (top panels) and selected ions 
(middle and bottom panels) in density-temperature phase 
space at z = 0.2 (left) and z = 1.0 (right). The grey pixe- 
lated shading is the fractional Q summed from SPH particles 
binned in 0.1 x 0.1 dex pixels. The temperature bimodality 
of metals described in £|3. 1 1 is obvious at both redshifts, with 
the metals pushing to lower densities and higher tempera- 
tures with time. 

We show the phase space occupied by Civ and Ovi 
absorption at z = 0.2 as the grey shading in the middle 
left and bottom left panels, respectively. These ions trace 
gas in a limited range in phase space, particularly weighted 
towards the cooler gas. Overlaid using coloured open circles 
are the ion mass density-weighted hydrogen densities and 
temperatures of simulated absorbers over a pathlength of 
Az — 14 (70 sight lines over z = 0.1 — 0.3). Colour and 
circle size both scale with the absorber column density as 
indicated by the legend above each panel. 

Our identified absorbers overlap with the darkest shad- 
ing, which provides a consistency check that our ion weight- 
ing accurately captures the physical conditions of the ab- 
sorbers. The two sets of histograms along the :r-axis (for den- 
sity) and y-axis (for temperature) directly compare the SPH- 
summed quantities (solid grey histograms) and the absorber- 
summed quantities (open red histograms); the amplitudes 
of these histograms can be directly compared. For example, 
the mismatch of C IV density histograms in denser halo gas 
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indicates that summing the halo C IV absorbers (red) un- 
derestimates the true ^Qj V arising from this phase (grey). 
The temperature histograms show that most of this miss- 
ing absorption is between io 40_41 K. The same is true for 
Si IV probing even higher densities (not shown). The main 
cause for this is that saturated absorbers tracing halo gas 
have their column densities underestimated without turbu- 
lent broadening, which is apparent in the differences between 
the vzw and vzw-turb CDDs in Figure [7] The mismatch is 
less severe but still noticeable for O VI. 

At z = 0.2, both C IV and O VI predominantly trace 
the cooler phases of gas at temperatures closer to 10 4 
K than 10 J K. Nv (not shown) traces the same tem- 
peratures, but at intermediate overdensities between C IV 
and O VI. For Ovi, this is in marked contrast to some 
previous work dCen fc Chisaril |201ll : ISmith et ail l201ll ; 
iTepper-Garcia et al.l 201 lh . allof whom find most of their 
O vi at T > 10 5 K at the same redshifts. We believe it is the 
interplay of our momentum-conserved winds with the cool- 
ing rates in the T ~ 10 5 K regime that results in our differing 
conclusion. We will show in i J6.2l th at using the pho to-ionised 
metal-line cooling rates of IWiersma et aL (|2009b) does not 
fundamentally change the phase space traced by Ovi ab- 
sorbers in our vzw model. Hence, w e reiterate the conclu- 
sion by lOppenheimer fc Pavel (|2009h that, for a physically 
well-constrained outflow model, Ovi primarily traces dif- 
fuse, photo-ionised metals. 

We show two higher ionisation species, Nevm and 
Si xil, at z = 1.0, because this is the redshift where 
COS is able to observe these lines. The Sixn doublet at 
499.3,521. lA falls into the COS sensitivity range at z > 1.1, 
but we show simulated absorbers and phase space centred 
at z = 1.0 because these characteristics are not strongly 
evolving. Sixn has a larger fraction of mass arising from 
collisional ionisation at densities corresponding mostly to 
halo gas. 

Surprisingly, most simulated Nevm absorbers are pre- 
dicted to be photo- ionised, tracing overdensities of around 
S ~ 10. However, the strongest simulated absorber is at 

10 5 ' 8 K gas, and 



N- 



Ne viii 



= 10 14 ' 2 cm~ 2 tracing T 



is likely more repre s entative of the three N e vill a bsorbers 
l|Savage et all 120051 ; iNaravanan et al l 120091 . l201ll ) already 
observed, since it has a similar column density and hence 
probably a similar environment. iMulchaev fc Chenl (|2009h 
find a 0.25 L* galaxy at an impact parameter of 73 /i -1 kpc 
at nearly the same redshift as the z — 0.207 Nevm ab- 
sorber along the HE0226-4110 sight line, suggesting that 
this absorption originates in collisionally ionised halo gas. 
They further suggest that the origin is a conductive interface 
between cool clouds moving through a hot medium, based 
on the alignment of Nevm with lower ionisation species, 
something that would n ot be captured in our simulations. 
INaravanan et all (|201lf ) note the existence of a 0.08 L* 
galaxy at 110 kpc at 5v = 180 krns" 1 from the Nevm 
absorber at z = 0.495 in the PKS0405-123 sight line, and 
further argue this is likely tracing T ~ 10 5 ' 7 K gas either in 
a hot halo or nearby WHIM. 

We emphasise that our model predicts numerous weaker 
Nevm absorbers tracing the diffuse photo-ionised compo- 
nent, as a separate and probably as yet undetected pop- 
ulation. These absorbers are mainly below -^Vj\f evm = 
lO 13 ' 7 cm- 2 (cf. 10 13 ' 85 , lO 13 ' 98 , lO 13 ' 96 cm" 2 for the dis- 



covered systems mentioned above) and should arise outside 
haloes. We cannot distinguish any other obvious observa- 
tional characteristics between the two types of Nevm ab- 
sorbers, as the simulated 6-parameters at COS resolution 
do not straightforwardly translate to temperature. Align- 
ment with Hi absorbers tracing overdensities of ~ 10, i.e. 
A^Hl ~ 10 13 ' 5 cm -2 at z = 0.75, is common for the photo- 
ionised Nevm absorbers. Upcoming deeper COS observa- 
tions might be able to detect this population. 

The true frequency of hot Ne vm absorbers may exceed 
our predictions given that our simulations do not adequately 
resolv e conductive interfaces within haloes. INaravanan et al.l 
(|2009r i estimates a higher frequency below z < 0.5, although 
their estimate is based on one absorber and is highly un- 
certain. Additionally, we predict a mass density of £^N e vm 
at z = 1.0 almost 3 times higher than at z = 0.2, which is 
suggestive of an even greater frequency of such absorbers at 
high redshift; the increase is primarily from the increased 
neon traced in the photo-ionised component. Given their 
distinct phase spaces (hot halo and dense WHIM vs. dif- 
fuse), the potential for absorber bimodality is greater than 
for O VI, but identifying it may require direct observation of 
an absorber's environment. We will return to this discussion 
in a future paper discussing the galaxy-absorber connection. 



5.1.1 Photo-ionised Absorption 

It may initially be surprising that high ionisation species 
such as Nevm, Mgx, and Sixn can be excited via photo- 
ionisation given their high ionisation potentials, but the ex- 
istence of such a bsorbers is straightforwa rdly traceable to 
the shape of the lHaardt fc MadavJ (|200lh background. We 
demonstrate this in Figure I12[ where we sho w the Line Ob- 
servability Index (LOX, iHellsten et al.l |T 998i )- the predicted 
equivalent width of metals lines associated with Hi if the 
IGM is uniformly painted with a metallicity of Z = 0.1 Zq 
and with a uniform temperature of T — 10 4 ' 2 K. Although 
these assumptions are overly simplistic, they are represen- 
tative and serve to illustrate the main point. 

At z = 1.0, Nevm peaks near nu ~ 10 _5 '° cm -3 
with equivalent widths as high as 8 mA. Our simulated 
photo- ionised Nevm absorbers are stronger than this ow- 
ing to the inhomogeneous enrichment, which results in 
higher metallicities for the detected absorbers, and further- 
more arise at slightly higher densities owing to the posi- 
tive inetalhcityj^densi^CT present in the vzw mod- 
els (|Qppenheimer fc Davdl2006l ). Photo-ionised Sixn should 
exist around the mean density (assuming the metallicity 
is sufficiently high), but it should not exceed even 1 mA, 
well below the detectability threshold of COS. Photo-ionised 
Mgx exists at higher overdensities but is still very weak. A 
detection of Mgx or particularly Si XII, while likely rare and 
exceedingly difficult, would provide the most unambiguous 
UV tracer of hot halo gas at T ~ 10 6 ~ 7 K. 



5.1.2 Ovi in Other Wind Models 

The phase space traced by Ovi at 2 = 0.2 is shown 
for the constant, slow, and no-wind models in Figure 1131 
While the cw model appears qualitatively similar to vzw 
with most absorbers being diffuse and photo-ionised, there 
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Figure 11. Q phase space diagrams at two redshifts (left- z = 0.2, right- z = 1.0) for all metals (top) and two common ion species 
observable with COS at each redshift (centre and bottom) from the r48n384vzw model. Grey shading indicates the fractional Q summed 
from all SPH particles binned in 0.1 X 0.1 dex pixels, and grey histograms show the distributions in density and temperature for the 
metals and each ion species. Overlaid coloured circles are simulated COS absorbers over Az = 0.2 for 70 sight lines where S/N=30. 
Colour and symbol size scale with absorber column density. Red histograms count the summed Qs of absorbers as functions of density 
and temperature, calibrated to the same scale as the grey histograms using Equation \3\ and are directly comparable to test how well 
simulated observations recover the total ion density contained in the simulations. Solid lines correspond to the division of the four phases 
as in Figure [2] and the dashed vertical lines indicates mean overdensity. 
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Figure 12. Line Observability Indexes (LOXs) at z = 0.2 and 1.0 
showing the photo-ionised equivalent widths for a variety of metal 
species as a function of H I column density assuming the density- 
Af Hl relation from Equatio n 3 of DIP, Z = . 1 Z w T = 10 4 ' 2 K 
for the metal species, and a lHaardt fc MadarJ j200ll') background 
renormalised to match the vzw Ly« forest flux decrement (D10). 



are significantly fewer absorbers at A^q vi > 10 14 ' 3 cm -2 . 
Stronger outflows create fewer strong absorbers, because 
metals preferentially reach lower overdensities where they 
cannot rapidly cool. Even though the IGM is generally hot- 
ter in this model relative to vzw, absorbers still remain 
predominantly photo-ionised. Despite the metal content of 
the cw WHIM being ~ 6.4 times higher than in the vzw 
model, the limited phase space traced by the O VI does not 
probe this enrichment, which can only be directly detected 
in species with X-ray transitions (e.g. Ovu and Oviii). 

Slow winds produce O VI statistics nearly indistinguish- 
able from the vzw model as reflected in the phase space 
plot. This same simulation produces a similar z ~ galac- 
tic stellar mass function to the vzw model below M* 
l|Oppenheimer et al.ll201Ch . thus the match with the Ovi 
statistics is consistent with sub-M* galaxies producing the 
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Figure 13. O VI Q, phase space diagrams as described in Figure fTTI 
for other wind models: constant winds (cw, f w i n d = 680 kms -1 , 
77 = 2), slow winds (sw, t'wind = 340 kms -1 , 77 = 2), and no 
winds (nw) for z = 0.2 O VI. 



Ovi absorbers (OD09). The no- wind simulation shows that 
winds are required to enrich the diffuse IGM; no absorbers 
exist in this phase. 

Overall, the conclusion that the O VI detected by COS 
(at least under the S/N assumptions here) will be predom- 
inantly photo-ionised is robust to variations in the wind 
model and the cooling rates (as we will show in i]6.2[) . Even 
Nevm, at columns just below what is currently accessi- 
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Table 2. Cosmic Densities of Ions 



Ion 


Redshift 


Model 


abs 


S2 SPH 


free 


Ovi 


- 0.5 


vzw 


36.7 


57.6 


64 




55 


cw 


20.0 


30.7 


65 


» 




sw 


41.6 


65.8 


63 


n 




11 w 


0.11 


0.30 


37 






vzw-turb 


45.9 


57.6 


80 


Civ 


- 0.5 


vzw 


11.5 


34.0 


34 






cw 


5.0 


12.2 


41 






sw 


8.5 


24.2 


35 






vzw-turb 


24.8 


34.0 


73 


Si iv 


- 0.5 


vzw 


2.3 


10.3 


22 






cw 


0.9 


2.6 


35 






sw 


1.3 


5.2 


25 






vzw-turb 


6.9 


10.3 


67 


Ne viii 


0.5 - 1.0 


vzw 


4.5 


10.7 


42 






cw 


4.0 


8.8 


45 






sw 


7.6 


15.5 


49 


55 




vzw-turb 


3.4 


10.7 


32 



a Summed ion cosmic density by SPH particle, in units of 10 — . 

6 Summed ion cosmic density by simulated absorbers, in units of 10" 

c Recovery percentage, i.e. n a b s /^sPH x 100. 



ble, becomes dominated by photo- ionised absorbers. Sixn 
at COS sensitivities probes hot gas, but not at low WHIM 
densities, and anyway these lines should be rare. A key con- 
clusion from this work is that the bulk of the diffuse collision- 
ally ionised WHIM gas, i.e. the so-called missing baryons, is 
only traceable through high ionisation X-ray transitions. We 
will explore such X-ray absorption in a forthcoming paper. 



5.2 Cosmic Densities of Ions 

The central theme in our exploration of the physical condi- 
tions of metal absorbers is the complex link between the en- 
richment patterns of metals and their observational tracers. 
Simulations show that even the most commonly observed 
low-z metal absorption line, O VI, traces only a small frac- 
tion of the metals a n d baryons outside of galax ies (OD09; 
ICen fc Chisarl l201ll ; iTepper-Garcia et alj 120111 ). Further- 
more, even when one uses Equation [3] to sum the observed 
absorbers to find the cosmic density of a particular ion, fi a bs, 
this does not necessarily yield the true physical mass density 
of that ion. 

We compare the mass-weighted cosmic ion density 
summed via the contribution of every SPH particle outside 
of galaxies, f2sPH (the grey shading in Figure to f2 a bs 
(Equation [3] the volume- weighted observational probe rep- 
resented by the coloured circles for individual absorbers). 
Ideally the two methods should yield the same answer; i.e. 
f2 a ba should recover the true Osph- We list the comparisons 
in Table [2] for a number of different ions and wind models 
along with the percentage of the SIsph recovered (/rec) in our 
simulated observational sample corresponding to S/N=30 
COS observations over a pathlength of Az — 35 and any 



species often have lower / rcc , which at least double for J7q iv 
and figjrv when turbulence is added to the vzw model. 
This is because these species (i) arise from higher densi- 
ties that are more stochastically sampled in our sight lines, 
and (ii) depend heavily on the derived column densities 
for the strongest absorbers, from where most of the cos- 
mic density of an ion arises. Siiv is the most stochastic 
absorber we follow, and arises from the highest densities 
(n H = 1CT 4 - 10~ 2 cm" 3 ). 

Conversely, £Iq vi and £l^e VIII are higher ionisation 
species that trace lower densities with less stochasticity. 
This explains the lower dispersion of recovery fractions be- 
tween the models. Recovery of ^Nevill (32-49% between 
z — 0.5 — 1.0) suffers because of sensitivity issues; the ma- 
jority of Nevm is not well-sampled at diffuse overdensities 
because it is too weak to be detected. Si XII is an even more 
extreme example with vzw having / reo = 13%; only a por- 
tion of the collisionally ionised Sixn halo component is re- 
covered. Ovi traces a statistically well-sampled phase space 
with the most sensitivity to provide the most complete cen- 
sus of f2sPH, but still manages to recover at most 80% of the 
, total. 

Another practical issue results from saturated absorbers 
that have uncertain column densities. AutoVP typically 
estimates lower column densities than the true value, and 
f2 a bs becomes under-estimated. Adding turbulent broaden- 
ing raises the column density where lines become saturated 
and, therefore, f2 a bs rises and also becomes more accurate 
as we discuss in i)6.3l This effect can be quite dramatic for 
species with high stochasticity such as C iv and Si iv at 
z < 0.5, where applying turbulent broadening more than 
doubles / reo . 

In summary, limited sensitivity and stochasticity can 
lead to inaccuracies in the determination of the cosmic 
ion abundance. Sensitivity issues always lead to an under- 
estimate, while stochasticity leads to more uncertainty in 
an ion's cosmic density. Additionally, saturated absorbers 
can hide significant amounts of true absorption, leading to 
under-estimates of cosmic ion densities, especially for low 
ionisation species. Hence while cosmic ion densities are a 
useful characterisation of the overall metal evolution in the 
IGM, detailed comparisons between models and data are 
more robustly accomplished using absorber statistics. 



6 MODEL VARIATIONS 

A fundamental difference between the Lya and metal-line 
forests is the latter's high sensitivity to outflow models. In 
this section we explore such sensitivities of the model pre- 
dictions to variations in some important physical and chem- 
ical assumptions. We will mainly focus on O VI, but will 
also consider other species. OD09 already explored a range 
of model variations for O VI, finding that this ion consis- 
tently traced the diffuse, T < 10 5 K IGM. G iven that other 

Cen fc Chisaril 



groups find different results for O VI (e.g. 



identified line below 10 15 cm 2 . This pathlength appears to l201ll ; ISmith et al1l201ll ; lTepper-Garcia et al.l 



Iml), 



adequately sample the densities from which the species we 
discuss arise. 

Between 22 and 80% of fisPH is recovered, depend- 
ing on the species and the wind model. Lower ionisation 



centrate our study on varying four aspects of the mod- 
elling: the ionisation background, metal-line cooling, turbu- 
lent broadening, and metal inhomogeneity. We present the 
EWDs and the CDDs for the variations that we explore in 



Metal-line absorption and physical conditions 21 



10 



50 100 



500 



13 



14 



15 



100 

A 

'n 10 

T3 



g 

T3 1 



100 - 

A 

^ 10 
T3 



> 
o 

G 

T3 1 



100 

A 

15" 10 

\ 
> 
o 

a 

T3 1 



100 

A 

10 

T3 



o 
G 

T3 1 



ovr z=o. 0-0.5' 




vzw 

vzw— a 1 .2 
vzw— a 1 .57 



vzw-a2.0 

Danforth (2008) " Vfv t 
Tripp (2008) 1> 4 



H 1 



(Wl z = 0. 0-0.5 J 



I I I I 



vzw 

■ vzw- 16 

vzw-PI 
cw-PI 




OVI z = 0.0-0.5 



vzw 



J^J/*--.-, vzw-turb 



cw-turb-PI 



H 1 1 J I I I I 



oil z = 0. 5-1.0 

vzw 

— ■ vzw-Osig 



H "t 1 



r: 



vzw-ZOl 



5 l> 



10 



50 100 
W (mA) 



500 



Figure 14. Equivalent width cumulative distribution plots for 
O VI systems at (2) = 0.25 for different model variations including 
the shape of the ionisation background (top panel, §6.1), differ- 
ent types of metal-line cooling (2nd, §6.2), turbulent broadening 
(3rd, §6.3), and the metal inhomogeneity (4th, §6.4). Models are 
described in the text. 



Figures [Til and [151 respectively, and provide cosmic Ovi ion 
densities at z < 0.5 in Table [3] 



o 




Tr 1 



o 



-1 



o 



O 



o 



O 



1 " OVI z = 0.0-0.5 



Danforth (2008) 



• vzw 

■ vzw— a 1 .2 
A vzw— a 1 .57 
Ovzw— a2.0 



OVI z = 0.0 



• vzw 
■ vzw-PI 

AcW-PI 




5—f- — 1 — 4- — 1 — 
.0-0.5 



• vzw 

Ovzw— turb 
Avzw-turb-PL 
Xcw-turb-PI 




• vzw 

* vzw-Osig 
x vzw-ZOl 




13 14 

log[N (cm" 2 )] 



Figure 15. Column density distribution plots for Ovi systems 
at (z) = 0.25 analogous to Figure [T4l 



6.1 Ionisation Background 

Our conclusion that photo-ionisation is responsible for the 
majority of the absorption by high-ionisation species is sen- 
sitive to the shape and intensity of the ionisation background 
at far-UV and soft X-ray energies. No direct observational 
constraints exist between th e observed slope of extreme UV 
(EUV ) spectra of quasars ijTelfer et all 120021: IScott et al.1 
l2004f ) and the X-ray background (jWorslev et al.ll2005t ). The 
ionisation potential required to photo-ionise Ovi at 8.4 Ry- 
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Table 3. Cosmic Densities of Ovi between z = — 0.5 



log[QovJ 



log[N ovl (cm 2 )] 



Model 


abs 


1J SPH 


/rcc 


vzw 


36.7 


57.6 


64 


vzw-c*2.0 


36.2 


53.4 


68 


vzw-ctl.57 


42.5 


75.3 


56 


vzw-cel.2 


41.5 


110.5 


38 


vzw- 16 


37.3 


56.6 


66 


vzw-PI 


26.2 


39.6 


66 


cw-16 


23.9 


35.6 


67 


cw-PI 


7.7 


11.5 


67 


vzw-turb 


45.9 


57.6 


SO 


vzw-PI-turb 


32.5 


39.6 


82 


cw-PI-turb 


11.5 


11.5 


100 


vzw-Osig 


34.4 


57.6 


60 


vzw-ZOl 


51.5 


71.4 


72 



a Summing ion cosmic density by SPH particle, in units of 10 . 

6 Summing ion cosmic density by simulated absorbers, in units of 10" 

c Recovery percentage. Ratio of f2 a ]-, s /f7sPH X 100. 



dbergs lies in this unconstrained gap (J. M. Shull, private 
communication) . 

We use the CUBA package |Haardt fc Madaull200j) to 
generate three alternative ionisation backgrounds varying 
the slope of the EUV portion of the quasar spectrum. Figure 
[T]displays these three other backgrounds, which are referred 
to as HM2005 backgrounds (after Haardt & Madau 2005) 
with a — 1.2, 1.57, and 2.0, where the spectrum above 1 
Ryd has an intensity F v oc u~ a . The HM2005 background is 
included as an input into the latest pack age of CLOUDY an d 
uses the a — 1.57 power law based on iTelfer et al.l (2002), 
but the range of slopes explored here is reasonable given the 
observed range and uncertainty in quasar spectral shapes 
(F. Haardt, private communication). The fiducial HM2001 
spectrum uses a = 1.8 and is most similar to the HM2005 
a = 2.0 spectrum. 

The hardest background, a — 1.2, has 16 times more 
intensity at the O VI ionisation potential relative to the soft- 
est background, a = 2.0, while leaving the background at 
the Lyman limit and, therefore, the Lycv forest nearly un- 
changed. Figure [TU] shows the phase space distribution of 
Ovi absorbers under these two extreme assumptions, and 
demonstrates that the peak density for £7q vi shifts from 
10" 5 1 to 10~ 4 ' 3 cm when moving from the softest to the 
hardest backgrounds (cf. grey a;-axis histograms). However, 
typical Ovi absorbers (Nq vi = lo 13 - 5 " 14 - 5 cm~ 2 ) do not re- 
flect such a dramatic shift. Recovery fractions of SIq vj also 
decline from 56% to 32% (Table [3]), mainly because Ovi ab- 
sorbers from haloes have underestimated column densities 
due to saturation. 

The top panel of Figure [15] displays the O vi CDD for 
the four backgrounds. As the EUV intensity is turned up, 
photo-ionised O vi moves to higher overdensities creating 
more strong absorbers and fewer weak absorbers; a power 
law fit to the CDD would give a shallower slope. The a = 1.2 
background model appears consistent with the EWD in 
the top panel of Figure [14] without the need for turbulent 
broadening, but the linewidths are still underestimated for 
strong absorbers compared to the observed 6-parameters 
(not shown). When a = 1.2, strong photo-ionised Ovi 
absorbers populate haloes while the a = 2.0 background 
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Figure 16. f! phase space diagrams as described in Figure 1111 
for z = 0.2 Ovi for the two extremes in our assumed ionisation 
background. The top panel shows the hardest background (a = 
1.2 for the EUV power law of quasars) and the softest background 
(a = 2.0). 



predicts barely any photo-ionised halo absorbers. If photo- 
ionised Ovi is clearly identified within haloes, it suggests 
that a hard meta-galactic ionisation background is pene- 
trating halo gas. 

Collisionally ionised O vi is a poor tracer of WHIM gas 
in all cases, because even the weakest ionising background 
efficiently photo-ionises oxygen in its collisionally ionised 
temperature range, peaking at T = 10 5 ' 45 K, to higher ion- 
isation states. This is why O vi fails to trace the copious 
amounts of metals pushed by cw winds into the low-density 
10 5-6 K WHIM; this oxygen ionises to Ovn and beyond. In 
summary, reasonable variations in the hardness of the EUV 
background are unlikely to alter our fundamental conclu- 
sion that most O VI absorbers detectable with COS arise in 
photo-ionised gas. 
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6.2 Metal-Line Cooling 

Even if the ionisation background photo-ionises most O VI 
at T > 10 5 K to higher states outside haloes, there still 
exists a potential for some O v i to trace WHIM tem- 
peratures. ICen fc Chisaril ^Olll ), ISmith et all (|201lh . and 
iTepper-Garcia et al l Hon]) all predict a significant compo- 
nent of collisionally ionised Ovi at overdensities as low as 
~ 10. One possible explanation is the efficiency of metal-line 
cooling in the ~ 10 K temperature regime. To explore this, 
we now compare two implementations of metal-line cooling- 
(i) the assumption of collisionally ionised equilibrium (CIE), 
and ( ii) cooling in the presence of the lHaardt fc Madaul 
|200lf ) photo- ionising background (PI). 

Our mai n simulations employ t he firs t implementation, 
namely that ISutherland fc Dopital (|l993l 'l CIE metal cool- 
ing rates are added to hyd rogen and helium rates calcu- 
lated in the presence of the lHaardt fc Madaul (|200lh back- 
ground. However, in the presence of a photo-ionising back- 
ground, the number of bo und electrons is redu c ed, le ad- 
ing to less efficient cooling. ITepper-Garcia et al.l l|201lf ) in 
thei r §5.3 demonstrate that photo-ionised cooling times us- 
ing IWiersma et all (|2009bl l rates in OWLS (Overwhelm- 
ingly Large Simulations), rates that include the effects of 
the UV background on the metal cooling, are significantly 
longer for 0.1 Zq and Z© gas at «h = 10 -5 cm -2 . This leads 
to a higher equilibrium temperature where cooling balances 
photo-heating. 

Photo-ionised rates are clearly the more physically cor- 
rect cooling rates, therefore we compare the results of us- 
ing our normal metal cooling implementation (adding CIE 
rates) to those using photo-ion ised metal cooling rates based 
on the|w lersma et al.1 <|2009bh solar abundance cooling ta- 
bles. To do so, we evolve four 16/i _1 Mpc simulation volumes 
at the same resolution as our 48/i _1 Mpc simulation volumes: 
vzw and cw winds with CIE cooling (vzw-16, cw-16) and PI 
(vzw-PI, cw-PI) cooling. 

We present the EWDs and CDDs for all but the cw-16 
model along with the 48 /i _1 Mpc vzw simulation in the sec- 
ond panels of Figures [Til and [T51 respectively. The vzw-16 
simulation is statistically identical to the larger vzw simu- 
lation for all O VI observations. This comes as no surprise 
because OD09 demonstrated that Ovi is mostly associated 
with low-mass, < 0.3M*, galaxies, which are well sampled 
in both volumes. This justifies using the smaller volume to 
explore these physical variations. 

Comparing the vzw-16 and vzw-PI simulation, there 
exists 70% as much O VI absorption in the latter, specifi- 
cally Qq vi , 3.26 x 10~ 7 versus 2.35 x 10~ 7 , with the reduc- 
tion evenly distributed in column density (Figure [151 middle 
panel). Examination of the Ovi phase space for vzw-PI in 
Figure [17] (upper left panel) reveals that most of the O VI 
absorption still arises from the diffuse O VI phase. However, 
there are subtle differences compared to the standard vzw 
model (bottom right panel). The locus of cool metals at 
8 > 10 (i.e. n H > 10" 5 ' 5 cm- 3 ) is no longer flat at ~ 10 4 K, 
below which CIE metal-line cooling no longer contributes, 
but declines from ~ 10 4 ' 5 K to below 10 4 K for increasing 
overdensity as seen also in OWLS (cf. [w iersma et al.|[201o[ 
their Figure 1). Ovi absorbers, therefore, trace somewhat 
hotter temperatures, but are still within the diffuse phase. 
Somewhat more O vi exists in the WHIM phase owing to 



this less efficient cooling, but still only a handful of identi- 
fied absorbers actually trace this phase. 

Examining the two different cooling models in the cw 
simulations reveals much larger differences (Figure [T7l right 
panels); f^Ovi declines by a factor of almost three using 
PI cooling, from 2.08 x 10" 7 to 7.1 x 10" 8 . The main dif- 
ference is that cw winds push metals to lower overdensities 
where photo-ionisation plays a larger role in retarding cool- 
ing. Less oxygen cools below 10 4 ' 5 K where photo-ionised 
Ovi ionisation fractions are the largest. 

As an aside, we note that galaxy growth remains rela- 
tively unaffected by this change. The amount of stars formed 
at z = is reduced by only 8% using PI cooling relative to 
CIE cooling. This mostly arises because the reduced metal- 
line cooling at lower densities hinders the recycling of w ind 
material back into galaxies l|Oppenheimer et alj 2010T ). A 
small reduct i on in overall star formation was also noted by 
I Smith et al.l (|201ll l when making the same comparison in 
their grid-based simulation (cf. the red and black curves in 
their Figure 1). 

In summary, including photo-ionised cooling rates does 
not alter our fundamental conclusions about the phases O vi 
and other absorbers trace in our favoured vzw model - O vi 
still traces mostly diffuse, photo-ionised gas. We also demon- 
strate t hat using a different wind m odel (cw), similar to that 
used bv lTepper-Garcia et al.l l |201lt ). produces a much larger 
change in O vi absorber properties when using photo- ionised 
metal cooling rates, because this model pushes many more 
metals into very diffuse gas where the photo-ionised cool- 
ing rates have a larger impact. We note that the vzw model 
does a much better job than the cw model in fitting the ob- 
served z = galactic stella r mass function below 5 x 10 10 Mq 
jOppenheimer et al.ll20l6T ). from which most O vi absorption 
arises (OD09). Hence, we favour the vzw model, although 
our tests do illustrate that the interpretation of observed 
O vi still could depend on the details of the input physics. 

6.3 Turbulent Broadening 

In OD09 we introduced a heuristic prescription to add sub- 
grid turbulent broadening to the ^-parameters of O VI, af- 
ter the conclusion of the simulation, using specexbin. This 
added turbulent broadening is a function of the physical 
density of the gas. The prescription adds an additional tur- 
bulent 6-parameter in quadrature where 

feturb = \A 4051 °g[ n H] 2 + 156741og[n H ] + 43610 kms" 1 , (4) 

over the range 71h = 10 -5 ' 31 to 10 -4 ' 5 cm -3 and 

foturb = 13.93 log(n H ) + 101.8 kms" 1 (5) 

over the range nu = 10 -4 - 5 to 10 -3 cm . Furthermore, 
&turb retains its maximum value of 60 kms -1 at higher densi- 
ties. See OD09 §4.4 for an in-depth discussion of the physical 
motivation for this model and its relation to the turbulence 
observed on kpc scales in C IV pa ired absorber s alon g lensed 
quasar sight lines at z ~ 2 - 3 l|Rauch et all 120011 ). OD09 
argue that the turbulence is related to outflows and decays 
on timescales of the order of a Hubble time. 

We already showed results for the vzw-turb model in 
Figures l6l through 1101 Here we combine it with the vzw and 
cw models with photo-ionised cooling rate (vzw-turb-PI, cw- 
turb-PI) shown in the third panels of Figures [14] and 1151 
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Figure 17. f2 phase space diagrams as described in Figure [TT1 for Ovi at z = 0.2 for simulations with photo-ionisational-dependent 
metal-line cooling (top panels) and the assumption of collisional ionisation for metal-line cooling (bottom panels). The left panels show 
the vzw models and the right panels show the cw models. We show the 16/i -1 Mpc volumes in all four cases, but the Ovi observables 
are statistically identical to the 48h~ 1 Mpc volumes for collisional ionisation cooling. 



Photo-ionised metal-line cooling rates reduce vzw O VI line 
frequencies and J1q vj by the same factor, one-third, either 
including turbulence or without it. The vzw-turb-PI model 
agrees well with the observed O VI CDD and EWD, and bet- 
ter than the vzw-turb model, which over-produces weak O VI 
absorbers. The vzw-turb-PI model is our favoured model for 
O VI as it has the most physically motivated combination of 
outflows, ionisation background, physical metal-line cooling, 
and turbulent broadening to explain both O VI and Lya for- 
est obser vations (DIP), as well as o t her observations of low- 2 
galaxies |Oppenheimer et al.ll2010l ; iDave et alJkoilal lbl). 

We show the C IV, Si IV, and Ne vm CDDs for this model 
along with the vzw-PI model in Figure [18] As before, our 
favoured model continues to overproduce the incidence of 
the lower-ionisation species. We do not believe this discrep- 
ancy can be mitigated by cosmic variance, because even our 
16/i _1 Mpc box is sufficient to sample the gas and galax- 
ies responsible for the absorption of these ions (cf. vzw and 
vzw-16 show nearly identical results). However, the idea of 
adding turbulent broadening for lower ionisation species that 
trace higher densities may not be as valid as in the case of 
Ovi. For instance, using lensed quasars probing the z ~ 2.3 
IGM. [Lopez et all (|2007l ) observed sub-kpc sub-structure in 



lower ionisation species including C IV and Si IV that did not 
exist in O VI. At low redshifts, all of the Si IV and most of 
the Civ in our models have densities where our turbulence 
model predicts 6 t urb > 40 kms -1 , but the observations of 
Civ suggest much narrower linewidths (e.g. ICooksev et al.l 
2010). Hence O VI may be tracing more diffusely distributed 
gas that is more affected by large-scale turbulence injected 
by outflows, while C IV and Si IV could be tracing individual 
high-density cloudlets. Regardless, models both with and 
without turbulence overpredict the C IV and Si IV absorp- 
tion. 

Figure fT^l shows the phase-space diagram for metals and 
absorbers in the vzw-turb-PI model. Si IV and stronger C IV 
absorbers trace much higher densities than O VI. This re- 
sults in many more strong C IV and Si IV components than 
in the vzw-PI model, because AutoVP finds stronger Ns 
for otherwise non-broadened saturated profiles and weaker 
absorbers are subsumed into the stronger ones. The his- 
tograms in Figure [19] show that these absorbers mostly arise 
in cool halo gas. A difficulty with our simulations is that 
we may not resolve coo ling instabilities in the halo gas (e.g. 
Mai ler fc Bullo ck 2004) that could result in much of the car- 
bon being in dense cold cloudlets, which would favour lower 
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Figure 18. Column density distribution plots for Civ, Si IV, and 
Nevill for vzw winds with photo- ionisation-dependent cooling 
rates and turbulent broadening. Error bars correspond to Poisson 



ionisation states than C IV. This is in principle easily possi- 
ble as only 3-5% of the carbon outside of galaxies resides in 
Civ between z = — 0.5 (vs. 13-15% at z — 2 — 4), and hence 
a small change in the ionisation rate at ~ 3 — 4 Ryd could 
dramatically change the C IV statistics. It will be interesting 
to compare our predictions of C II and C ill to observations 
when a sufficient sample becomes available, to see if we sig- 
nificantly underpredict the amount of carbon in these states 

corresponding to our overpredicti on of C IV. 

Another possibility is that the lChieffi fc Limongil (|2004l ) 
Type II carbon yields we use are too high. For example, the 
IWooslev k. Weaver! (|l995T l carbon yields are ~ 0.2 dex lower, 
which could explain some of the differences with respect 
to the observations. The AGB stellar yields we use cannot 
explain the over-estimate, because carbon production from 
AGB st ars is sub-dominant to Typ e II SNe in our simu- 
lations (|Qppenheimer fc DaveHioog ). Yet another potential 



explanation is that a substantial fraction of carbon and sili- 
con are depleted ont o dust, while the more volatile oxygen is 
not. IZu et a.1.1 (|201lT) find that the vzw model c an explain the 
obser vational evidence for intergalactic dust l| Menard et al.1 
I2010T ) if the IGM dust-to-gas ratio is about half that in the 
Milky Way ISM, which about 25% of intergalactic metals 
(by mass) in dust. 

Ne viii traces similar overdensities as O VI and hence 
should be affected by the same turbulence, and hence the 
vzw-turb-PI model should be the most appropriate. But as 
Figure [18] shows, the predicted column density distributions 
do not vary significantly with PI cooling or turbulence. In 
detail, between z = 0.5 — 1.0, ^N e viii slightly increases with 
PI cooling, but declines by a factor of 20% with the addition 
of turbulence to 4.0 x 10~ 8 , because already weak lines are 
broadened enough so that AutoVP fails to identify them. 
Most of the Nevm absorption still traces photo-ionised gas, 
but the reduced cooling raises the bulk of Nevm absorbers 
from T = 10 4 " 4 ' 4 K to lo 4 - 4 " 4 9 K and leads to 22% com- 
pared to 7% of the Nevm absorption arising from T > 10 5 
K gas (Figure [TUJ) . This percentage could grow even larger if 
our simulations adequately resolved interfaces between cool 
and hot gas within haloes, which they do not as we discussed 
in SgU 

We show the cw-turb-PI model to demonstrate that tur- 
bulent broadening affects other wind models in much the 
same way as the vzw model, by making more wide O VI ab- 
sorbers and slightly fewer weak absorbers (cf. cw-PI model). 
This more appropriate cooling rates applied to the cw model, 
combined with turbulent broadening to match observed line 
widths indicates this model far under-produces Ovi, espe- 
cially weaker systems, when compared to observed EWDs 
and CDDs. 



6.4 Metal Inhomogeneity 

The degree of metal inhomogeneity in the IGM remains 
highly uncertain, but dynamically very important. The 
distribution of metals affects gas cooling, which is criti- 
cal for galaxy formation because the temperature of the 
gas determines the e f ficiency of gas accretion onto galax- 
ies jKatz et al. 2003; iKeres et al.1 |2005| ; iDekel fc Birnboiml 
120061 ; IKeres et al.l 12009 ). More metal mixing leads to more 
star formation as more efficient met al-line cooling affects 
more baryons (Wie rsma et al.l l2009af ). In our simulations, 
we do not mix metals, and each wind particle retains its 
own metals as it propagates out of galaxies. 

The homogeneity of metals impacts the statistics of ab- 
sorbers, which we explore here. We consider two ad hoc 
models for distributing metals throughout the IGM, by re- 
placing particle metallicities in our simulation with prescrip- 
tive metallicities in our vzw model. In the first (vzw-ZOl), 
we assume a uniform metallicity of 0.096 Zq in the IGM 
(vzw-ZOl), which is the global z = IGM metallicity in 
our vzw model. Note that this is similar to 0.1 Zq painted 
onto earlier simulatio ns to reproduce older observed Ovi 
absorber frequen cies l|Cen et al.1 l200ll ; iFang fc Brvanl 1200 ll ; 
IChen et al.ll2003T ). The second model ( vzw-Osig) takes the 
average [O/H] of the vzw model in each 0.1 x 0.1 dex pixel 
of phase space (upper left panel of Figure [20] at z = 0.25) 
and paints this metallicity onto all SPH particles within that 
phase space cell, with zero dispersion. Both prescriptions 
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Figure 19. f! phase space diagrams of the vzw-turb-PI model for metals, Si IV, Civ, and Ovi at z = 0.2; plus metals and Nevill at 

2=1. 
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preserve the total metal mass, but neither model is dynam- 
ically self-consistent, because metal cooling would alter the 
physical state of the gas and feedback would have to enrich 
the baryons in the first place. However, our purpose here 
is to explore extreme cases of metal homogeneity to see if 
absorption line statistics can distinguish the level of metal 
inhomogeneity. 

Figures [14] and 1151 bottom panels, show the EWDs 
and the CDDs of these models, compared to our original 
vzw model. The vzw-ZOl model, as found in OD09, pro- 
duces too many weak absorbers, too few strong absorbers 
(in terms of column density), and a large alignment fraction 
with Hi (Figure [21) . Collisionally ionised Ovi (T > 10 5 K) 
accounts for a quarter of Hq vi and two-thirds of lines above 
Nq vi > 10 14 cm -2 . Simply painting metals onto a simula- 
tion uniformly overestimates O VI and other species in CIE, 
because the most commonly observed UV resonance lines 
(Civ, Ovi, and soon Nevm) are also some of the most im- 
portant coolants. 

Painting the average metallicity by phase space pixels 
explores how the dispersion in metallicity among SPH par- 
ticles affects the line statistics. This model shows negligible 
differences versus the original vzw case, showing that the 
fact that metals are not mixed in our simulations has no im- 
pact on Ovi absorber statistics. We reiterate that this test 
does not dynamically account for differences in the cooling 
rates that would result from such metal mixing, and hence it 
is not self-consistent; this test isolates the impact of spatial 
inhomogeneities in the metals. 

Figure [201 provides a broader view of metal inhomo- 
geneity in our original vzw simulation. The upper left shows 
the oxygen metallicity in phase space, this is similar to Fig- 
ure [2] except here we only consider oxygen and the general 
trends are as before. The upper right panel shows the frac- 
tion of baryons enriched with oxygen, / cn , using a threshold 
of [O/H] = 10~ 6 Zq. / on is not very sensitive to this thresh- 
old and is used to discard SPH particles enriched to very 
low levels by delayed feedback. Practically all the ICM gas 
and most gas, both hot and cold, inside of haloes is enriched. 
Interestingly, the / en distribution for the nw simulation (not 
shown) looks nearly identical to the vzw model above 10 s K, 
indicating that enrichment mechanisms other than outflows 
(e.g. dynamical stripping, delayed feedback) are getting hot 
baryons, i.e. hot halo and WHIM gas, above this threshold 
metallicity. In the IGM outside of haloes, most particles re- 
main pristine, as the typical fractions of enriched particles 
are < 0.1. In photo-ionised gas within voids, the fraction of 
enriched particles can drop below 10 -3 . 

The lower left panel shows the average metallicity of 
only the enriched particles, [0/H] en . This demonstrates 
that much of the enriched portion of the diffuse IGM has 
[O/H]=0.1 — 0.3 Zq. Since the outflows carry out metals 
that are representative of the ejecting galaxy, this indicates 
that the dominant galaxies enriching the IGM tend to have 
sub-solar metalliciti es; these galaxies are genera lly well be- 
low L* (as found in lOppenheimer fc Pavel lioosh . However, 
there is a selection effect that highly enriched IGM parti- 
cles undergo stronger metal-line cooling, dropping them out 
of the warmer phases onto the photo-ionised locus. Hence 
enriched particles in the underdense regions actually have 
higher metallicities than the WHIM gas in moderately over- 
dense regions, since these are the only particles that can cool 



at such low densities. Nevertheless, the / cn shows that such 
particles are quite rare. 

The lower right panel shows the log-normal dispersion, 
(j([0/H] cn ) of the metals around t he median log[ | Q/Hle n1 for 
the enriched partic les only. Both ISchave et alj (2003) and 
ISimcoe et al.l (|2004[ ) found that C iv and O vi enrichment at 
z ~ 2 — 4 could be described by a log-normal distribution, 
but applied to all baryons. In general, regions with disper- 
sions below 0.5ct are enriched via outflows consisting of ISM 
gas having similar metallicities. The higher dispersions cor- 
respond to regions enriched via other mechanisms such as 
dynamical stripping and delayed feedback from intracluster 
stars. 

6.5 Alignment Statistics 

The association of H I with O vi offers a potentially powerful 
way to explore metal inhomogeneity. We show alignment 
statistics of Ovi components with nearby Hi components 
in Figure [21] We plot the fraction of Ovi absorbers with an 
associated Hi absorber at or below the velocity difference 
indicated on the s-axis, with the criteria that both absorbers 
have W ^ 30 mA and z < 0.5. Overlaid are STIS data from 
iTripp et alj (|2008l ) using the same selection criterion. 

The data are unable to discriminate between vzw and 
vzw-Osig, but disfavour the vzw-ZOl model, which tends to 
produce more alignment. OD09 showed that Ovi absorbers 
trace similar diffuse densities and temperatures as Hi, but 
that the Ovi-traced enriched gas is inhomogeneous and 
clumpier than the diffuse baryons traced by the Lya forest. 
The difference in metal inhomogeneity between the vzw and 
vzw-Osig models does not show up in the EWD and CDD 
distributions in Figures [141 and [TBI but the expected obser- 
vational sample of the two most common species probed by 
COS is expected to provide the ability to statistically dis- 
tinguish these two models. The cw model shows lower align- 
ment owing to faster winds pushing metals further from the 
Lya in both physical and velocity space. Both turbulence 
and photo-ionisation-dependent cooling (not shown) do not 
alter the statistics from the vzw model. 

Alignment statistics versus H I fundamentally probe the 
level of photo-ionisation in metal absorption, because H I 
arises almost exclusively in photo-ionised gas. Models that 
produce more collisionally ionised O vi will therefore have 
significantly less alignm ent with H I. Curre nt observations of 
Ovi-Hl alignment from lTripp et alj (120081 ) shown in Figure 
[2TJ suggest that the vzw model produces roughly the correct 
level of photo-ionised gas in relation to Hi, and therefore 
strongly favours our conjecture that much of the Ovi seen 
with COS is photo-ionised. 



7 DISCUSSION 

Our favoured model includes momentum-conserved winds 
and uses photo-ionisation-dependent cooling rates. The vzw- 
turb-PI model matches the Ovi statistics including the 
z < 0.5 EWD, CDD, fiQ VI , and 6-parameters when tur- 
bulent broadening is added as argued and theoretically mo- 
tivated in OD09. We further advocate vzw winds for their 
ability to provide an adequate fit to th e z = galactic stel- 
lar mass function below 5 x 10 10 Mq l|Qppenheimer et al.l 
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Figure 20. Phase space diagrams displaying oxygen metallicity characteristics for the r48n348vzw simulation at Z = 0.25. The colour 
scales correspond to mean metallicity (upper left), the fraction of baryons enriched to at least 10 -6 Zq (upper right), the median 
metallicity of the enriched fraction (bottom left), and the log- normal metallicity dispersion of the enriched fraction (bottom right). 
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Figure 21. Alignment fractions for Ovi with Hi components 
between z = — 0.5 as a function of inter-species component 
velocity difference. Both Ovi 1032A and Hi 1216A are required 
to have W > 30 mA. 



l2010t ). from which most Ovi absorption arises (OD09), as 
well as other properties o f the low-redshift galaxy pop- 
ulation (|Dave et al.l l2011al fbT); these galaxy properties are 



not expected to be greatly impacted by adding IGM 
turbulence or PI cooling. The HM2001 ionisation back- 
ground multiplied by 1.5 appears capable of both re- 
producing the observed Lycv forest (D10) and Ovi that 
is principally photo-ionised. The vzw wind model also 
has successes fitting higher redsh ift data, including th e 
z — 2 mass- metallicity relationship (IFinlator fc Daveli2 008). 
the z ^ 6 galaxy luminosity function (|Dave et alj 1200a ; 
IFinlator et al] I2010T ). and observations of IGM metal en- 
richment at redshifts ^ 1.5 dOppenheimer fe Davi 1 20061 : 
lOppenheimer. Dave, fc Finlatorl l2009f ). The most obvious 
failing of the vzw model is its inability to quench star forma- 
tion in > M* galaxies and create the observed red sequence, 
although this likely has to do with our simulations lacking 
a mas s-dependent quenching m echanism such as AGN feed- 
back l|Gabor et alj |20ld . l201lh . However, massive galaxies 
are not responsible for the vast majority of O VI enrichment 
(OD09), because they enrich a comparatively small volume 
of the Universe, and so this quenching of massive galaxies 
is not expected to significantly impact our IGM enrichment 
levels. 

Nevertheless, this model also has significant problems 
that are shared by the other wind models we consider here. 
C IV and Si IV observations pose a problem in that our 
favoured model produces t oo many absorbers compared to 
ICooksev et ail (|2010|J201 il l . It might be possible to alter the 
ionisation background at 3-4 Rydbergs to mitigate this dis- 
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crepancy while still satisfying the other observational con- 
straints. Updated ioni sation backgrounds gene rated with the 
new version of CUBA (|Haardt fc Madaull201ll ) are worth ex- 
ploring in the future. Alternatively, our assumed yields could 
be too high, or car bon and silicon could be preferentially 
depleted onto dust (jZu et al.ll201ll ). However, the discrep- 
ancy may result from limitations present in our simulations, 
namely that we fail to resolve cooling instabilities in halo 
gas where such lower-ionisation absorbers are expected to 
arise. Quantifying C n and C ill absorption may help distin- 
guish between these possibilities. Our favoured model also 
under- predicts the frequency of Nevill li nes with column 
densi ties similar to those observed so far l|Naravanan et all 
2009), which again may be a result of insufficient resolu- 
tion within the halo gas, in the sense that we fail to resolve 
metal-bearing interfaces between cold clumps embedded in 
hot, tenuous halo gas. Relatedly, OD09 (§5.4) found a lower 
frequency of complex, collisionally ionised O VI absorbers, 
which may also arise in shock interfaces within hot halo gas, 
as implied by some observations. 

In contrast, absorbers arising in more diffuse gas should 
be more robustly modelled in our simulations. Hence, we 
are confident that, when probed to sufficiently low column 
densities, most high ionisation lines of Ovi and Nevill will 
mostly trace diffuse photo-ionised gas. Since there have been 
a number of recent publications exploring Ovi and other 
metal-line statistics in cosmological simulations to z = 0, 
we now compare and contrast our results with these works. 

7.1 Feedback Methods in Other Simulations 

Our discussion of alte rnative feedback methods begins with 
ICen fc Chisaril (|201ll ). who generate statistics for O VI and 
Civ from 2 = 6 -> in a 50 /i _1 Mpc 2048 3 -cell uniform 
mesh simulation. Their winds are implemented b y distribut- 
ing g as and energy onto the 27 nearest cells l|Cen et al.l 
120051 ). which instantaneously distributes metals over a 73 
comoving /i -1 kpc cube. Clearly this model does not self- 
consistently track the dynamics of enriching this cube, but 
it nevertheless provides a heuristic way to distribute met- 
als into the IGM. Their favoured wind efficiency is egsw = 
7 x 10~ 6 , where eGsw represents the fractional rest mass 
energy of all stars formed put into thermal feedback. Their 
total wind energy is t hen -E w i n d = 1-2 x 10 49 erg/ Mg, which 
ICen fc Chisarl (20 111 ) argue is roughly comparable to the to- 
tal SNe energy in a Ch abrier initial mass function (IMF). 

ISmith et all l|201lh similarly use fixed-mesh Eulerian 
simulations and parametrise star formation and feedback the 
same way. They can reproduce much of the observed z < 0.4 
O VI CDD with a wind energy efficiency of egsw = 10 -5 . We 
calculate their total -E w i n d = 1.8 x 10 49 erg/ Mq, which is 
nearly twice the SNe energy available in a Chabrier IMF if 
all stars 10 Mq explode as 10 erg S Ne, but about equiv- 
alent if all stars ^ 6 Mq explode as SNe. ISmith et all (|201ll ) 
ad vocate the need of di stributing feedback over 27 cells (as 
m ICen fc Chisaril l201ll) to fit the observed star formation 
rate density history, whereas limiting the feedback to one 
cell makes it much less efficient at curtailing star formation. 
This means that in their case, metals are instantaneously 
distributed over a 99 to 294 /i _1 kpc cube depending on their 
simulations. 

Our simulations, in contrast, attempt to directly mimic 



observations of galaxies that indicate wind s peeds of hun- 
dreds of km/s emanating from galaxies (e.g. ISteidel et al. 
2010l). and scaling with the circular velocity (e.g. iMartin 
2005al ). We are inevitably limited by resolution and numer- 
ics in that we cannot fully resolve how winds escape from 
galaxies, so we are forced to employ parametrisations for 
these winds. The adopted parametrisations (including hy- 
drodynamic decoupling) are intended to represent observed 
winds at the scales at which they are observed. It is inter- 
esting that our winds typic ally travel ~ 100 physical kpc 
|Oppenheimer fc Dave 2008), which is similar to the sc ales 
over which lCen fc Chisaril (|201ll ) and lSmith et all (|201ll ) as- 
sume outflows distribute gas. The difference is that some of 
our winds can travel much farther or not as far, depending on 
the self-consistently tracked interaction between the gas and 
the ambient medium; we do not directly parametrise how far 
winds and metals travel. Moreover, in these fixed mesh sim- 
ulations the distance metals travel is explicitly governed by 
the spatial resolution, while in our case the wind speeds are 
based on the galaxy mass, which is a well-converged quan- 
tity. For instance, the lower-resolution mesh simulations of 
ICen fc Fan e (2006) produce a significant Ovi signal with- 
out any mechanism for distributing gas to nearby cells, in 
contrast to our no-wind case that produces essentially no 
IGM metal absorbers. The direct sensitivity of fixed mesh 
simulations to ad hoc numerical parameters makes it diffi- 
cult to interpret their conclusions about the physics driving 
IGM enrichment. Adaptive mesh refinement (AMR) simu- 
lations could potentially provide a less resolution-sensitive 
implementation of winds. 

Finally, we point out that our feedback model has quite 
modest energ y requirements relat ive t o the above mesh sim - 
ulations. Both lSmith et al.l l|201ll) and lCen fc Chisaril (|201ll ) 
strain the limits of available sup ernova energy. Our cw mode l 
and the feedback prescription in lTepper-Garcia et alj (|201ll ) 
("Uwind = 600 kms -1 , 77 = 2) require substantial fractions 
of SN energy to be converted directly to kinetic winds- 
Ewnd/EsN ~ 0.75 — 0.97 assuming the (conservative) 10 Mq 
lower mass limit for SNfl Comparatively, the vzw model 
has -Ewind/-EsN = 0.56 at z = 1 and 0.45 at z = 0, which 
is meant as an illustrative comparison since this wind pre- 
scription simulates energy generated during the lifetimes of 
O and B stars(3 Because vzw feedback energy scales with 
a, the lower mass galaxies that are primarily responsible 
for the enrichment of the diffuse IGM have even less of an 
energy requirement. For instance, a M* = 10 10 Mq galaxy 
has g W md/£sN = 0-30 at z = 1 and 0.21 at z = 0. Note 
that IStrickland fc Heckmanl (|2009l ) in their §5.4 misinter- 
pret the vzw model as requiring more energy than in the 
SN energy budget; this is true only for central galaxies in 
haloes more massive than 10 13 ' 5 Mq at z — 1, or for between 
2-4% of all wind particles launched at any redshift. These 
galaxies are not responsible for enriching the IGM, and in 

3 The iTepper-Garcia et al. 1120111) model converts to 
^wind/^SN = 0.40 using the 6 Mq lower mass limit for 
SNe, and the cw model gives an efficiency of 0.52 using this 
limit; however we scale all models to the 10 Mq mass limit in 
our comparison. 

4 The E wind /E SN = . 72 r48 n384vzw value published in Table 1 
of lOpp cnhcimer et al.l ll2010t ) is a miscalculation, overestimating 
the efficiency by 50%. 



30 B. D. Oppenheimer et al. 



any ca se should not have significant ongoing star formation 
today (|Gabor et alj|20ich . 



7.2 Metal Mixing 

Absent in our simulations is any attempt to simulate the 
mixing of metals. For the purposes of cooling and generating 
mock s pectra, we use particle meta llicities, which is in con- 
trast to iTepper-Garcia et alj (|201ll ). who smooth metallici- 
ties over SPH particle smoothing lengths when calculating 
metal-line cooling. The ir metal smoothing was introduced in 
W iersma et alj (|2009al ) , who calculate that smoothed metal- 
licities are ~ 2/9th the individual particle metallicities when 
all surrounding particles in their kernel are pristine. Given 
that our metal-enriched IGM wind particles are often sur- 
rounded by pristine gas, the cooling times of our metals 
could be as much as ~ 4.5 times shorter using our indi- 
vidual particle metallicities if metal-line cooling dominates 
total cooling. We suggest this could be the source of the 
greate st differences with the results of ITepper-Garcia et al.l 
l |20 111) , who find the majority of O VI above 10 5 K. Our cw- 
PI model is most similar to their "REFERENCE" model 
with notable other differences including a slightly higher 
wind speed (v W md = 680 vs. 600 krns" 1 ) and that we decou- 
ple our winds; both of these should make it easier for winds 
to escape galaxies in our simulations. Nonetheless, both sim- 
ulations have similar Qq vi and similar overdensities bearing 
Ovi gas, but most of our Ovi gas resides around 10 4 ' 5 K, 
whereas theirs is at io 4 - 8-5 - 5 K with significantly more Ovi 
at temperatures > 10 5 ' 1 K. 

Fixed mesh simulations with distributed feedback also 
can lead to much smoother metallicities, resulting in longer 
cooling times and hotter metal-bearing gas. We argue that 
the lack of metal smoothing is a fundamental difference caus- 
ing a significant photo-io nised component of Ovi in our 
simul ations versus others (ICen fe Chisarj 2011 1_ Smith et al 



l201ll;lTepper-Garcia et al.H201ll ). Both lTepper-Garcia eta! 
1 201ll ) and lSmith et al.l (|20 1 if ) make a strong case that us- 
ing CIE metal-line cooling rates over-estimates cooling re- 
sulting in too much photo-ionised Ovi in OD09. This is no 
doubt true, especially given the fact that in our cw simu- 
lations, photo-ionised metal-line cooling rates reduce Qq vi 
by a factor of three. But as we demonstrated in H6.2I for our 
favoured vzw case the difference is much less. Still, metal 
mixing can be an important determinant for the IGM phase 
of O vi absorbers. 

Treatments of met al diffusion have been implemente d 
into other simulations (|Greif et al.ll2009l ; IShen et alj|20ld) . 
but it remains unclear to what extent metals mix in the 
IGM. Interpreting observations of metal absorbers at higher 
redshift suggest that metals could remain in concentrated, 
smal l cloudlets w i th metallicities similar t o galactic ISM 
gas (|Simcoe et al.l 120061 : fechave et ai]|2007l ). This scenario 
would favour our lack of mixing and the fact that winds 
carry ISM gas directly into haloes and the IGM. In sum- 
mary, the necessity and specific implementation of metal 
mixing remains uncertain when considering the low-z IGM. 
It appears that metal mixing is as determinative as imple- 
mentations of metal-line cooling or feedback for interpreting 
IGM metal-line observations. 



8 CONCLUSIONS 

We examine the absorption and physical conditions of met- 
als residing in the ~90% of baryons outside of galaxies 
between z = and 2. We employ 48/i _1 Mpc, 2 x 384 3 - 
pa rticle cosmological simulatio ns, the same simulations used 
by lOppenheimer et al.l (|20 10T ) to explore the 2 = galac- 
tic stellar mass function and by iDave et alj (|201p| ) to ex- 
plore t he Lyq forest. Our work bu ilds on the previous re- 
sults of lOppenheimer fc D ave (2009), finding fundamentally 
the same conclusions for O VI, and also makes predictions 
for Si iv, Civ, Nevm, Mgx, and Sixn specifically tailored 
to the expected instrumental characteristics of the Cosmic 
Origins Spectrograph (COS). We subdivide the gas phases 
between the true IGM (5 < 120 at z — 0) and gas within 
dark matter haloes (120 < 8 < 10 5 ' 7 at z — 0), and between 
hot (T > 10 5 K, WHIM and hot halo) and cool (T < 10 5 K, 
diffuse and condensed) phases. 

Metal-line absorption is fundamentally different than 
its primordial big brother, the Lya forest, owing to its much 
greater sensitivity to the processes of galaxy formation and 
evolution. While galactic outflows suppress star formation 
and therefore metal production, they concurrently distribute 
metals into haloes and the IGM. Our simulation without 
winds (nw) produces 2.1 times more metals by z — than 
the momentum-conserved wind (vzw) model, but the vzw 
simulation has 3.5 times more gaseous metals outside of 
galaxies, and 380 times more metals in the diffuse IGM 
phase from which most IGM metal absorbers arise. The re- 
sulting enrichment patterns are quite sensitive to the details 
of the enriching outflows. The combination of wind velocities 
and environment determines to how low a density the winds 
reach, e.g., whether or not metals leave their dark matter 
haloes. Low velocities from dwarf galaxies easily reach the 
IGM, while > 1000 kms -1 winds from > M* galaxies of- 
ten fail to escape their haloes because of the greater am- 
bient density of the surrounding gas l|Oppenheimer et al.l 
2010). The wind velocity also determines the temperature 
to which the gas shock heats as the wind decelerates. Fast 
winds that shock gas to > 10 6 K are more likely to join the 
WHIM or hot halo gas. The density and temperature where 
the metals are deposited in turn impacts their cooling rate, 
which can return enriched gas to photo-ionisation tempera- 
tures. Tracking all these processes dynamically within grow- 
ing large-scale structure presents a complex and challenging 
problem that impacts our understanding of a wide range of 
phenomena related to the cosmic evolution of baryons. 

Our main results are as follows: 

• Metals progressively are found in more overdense re- 
gions from z = 2 — >■ 0. By z = in our favoured momentum- 
conserved wind (vzw) model, only 4% of metals are out- 
side of haloes (i.e. at densities below the threshold de- 
fined by Equation [1} , compared with 20% at z = 2. The 
constant wind (cw) model has many more metals outside 
of haloes, but they mostly reside in the difficult-to-detect 
WHIM phase. The majority of metals today reside within 
galaxies, in contrast to the bulk of baryons that reside out- 
side of haloes (§3). 

• For all the wind models, the mean metallicity of most 
of the baryonic phases is very slowly rising since z ~ 2. The 
exception is the condensed phase, which shows a more rapid 
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rise as winds progressively enrich halo gas closer to galaxies 
at later epochs. By z = 0, the IGM metallicities in both the 
cool (diffuse) and the hot (WHIM) phases are ~ 1/50 Z© 
in the vzw model, and Z = 0.096 Zq in the 90% of baryons 
outside galaxies (§3). 

• The IGM is enriched in an "outside-in" fashion where 
lower overdensities are enriched at higher redshift. Higher 
ionisation metal species generally trace older metals at lower 
overdensities owing to these species being primarily photo- 
ionised. The w w i n< j oc a relation for vzw outflows results in 
a continuously rising metallicity-overdensity gradient and 
a gradual age-overdensity anti-correlation, which contrasts 
with the constant wind simulations (§3.3). 

• Constant (cw), momentum-conserved (vzw), and slow 
(sw) wind models all predict z < 0.5 Ovi observational 
statistics (equivalent width distributions, column density 
distribution, cosmic ion densities) that agree within a fac- 
tor of two with observations by STIS and FUSE. By con- 
trast, the no- wind model produces an Ovi incidence 100 
times lower than the other feedback models and observa- 
tions, which strongly emphasises the need for galactic su- 
perwinds to enrich the IGM to the observed levels. All wind 
models over-produce the C IV and Si IV incidence relative to 
available observations (§4). 

• O VI, C IV, and Si IV absorbers predominantly (by 
number) arise in photo-ionised, cool gas. Photo-ionisation 
results in higher ionisation lines tracing lower gas densities: 
Si IV and C IV arise primarily from halo gas at z < 0.5, while 
Ovi arises mainly in the diffuse IGM (§4,5.1). We predict a 
population of Nevill absorbers with ^VjvJevill ~ 10 13 ' 6 cm~ 2 
tracing photo-ionised diffuse T < 10° K gas at overdensities 
of ~ 10. Stronger Nevill absorbers (like stronger Ovi sys- 
tems) are more likely to be associated with T = io 5 ' 5-6 
K gas and may share more in common with the absorbers 
already detected (§4,5.1). These results do not change when 
using metal-line cooling rates including the presence of a 
photo- ionising background (§6.2,6.3). 

• Mg x and Si XII absorbers are very rare, with frequen- 
cies of dn/dz « 1 for COS S/N=30 spectra in all of our 
models. These absorbers likely trace halo gas at T > 10 6 
K (§4,5.1). These species and Nevill all have a significant 
photo -ionised componen t , whic h is consistent with the shape 
of the lHaardt fc Madaul (|200lh ionisation background; how- 
ever Mg x and Si XII are too weak in this state to be detected 
using COS (§5.1.1). No UV tracer we examine consistently 
and predominantly traces WHIM gas, i.e. shock-heated gas 
outside of galaxy haloes. 

• Using metal-line cooling rates accou nting for photo- 
ionisat ion by the ionisation background (Wi ersma et al.l 
2009b) reduces the incidence of O VI absorbers at all column 
densities by one-third in the vzw model, which still provides 
an acceptable match with the observations. The difference is 
more dramatic for cw winds since it pushes metals to lower 
overdensities; Ovi statistics are reduc ed by two-thirds and 
arise from gas closer to T = 10 s K. iTepper-Garcia et al.l 
l|201lf) and ISmith et all ( |201lf > correctly argue that assum- 
ing collisional ionisation equilibrium for metal cooling rates 
(with no photo-ionisation) leads to over-efficient metal-line 
cooling; however we find that the moderate vzw winds leave 
metals at higher overdensities where metal-line cooling is 
still efficient even in the presence of photo-ionisation and 



results in most of the O VI absorbers tracing T < 10 5 K gas 
(§6.2). 

• We reiterate the result of OD09 that heuristic tur- 
bulent broadening is required to reproduce the observed 
Ovi linewidths, but we suggest that this prescription may 
not be applicable to C IV and Si IV primarily arising from 
haloes. O VI arises from longer columns with potentially 
many cloudlets with random/turbulent velocities at sub- 
particle scales, while lower ionisation species are more likely 
to trace individual higher density cloudlets (§6.3). 

• Absorber alignment statistics of H I with O VI provide 
an independent constraint on the level of metal inhomogene- 
ity in the IGM, which we quantify in §6.4. We disfavour a 
uniform metallicity of Z = 0.1 Zq from comparisons with 
existing data, and future COS observations can potentially 
quantify the fraction of enriched gas and the dispersion of 
metals in the IGM (§6.5). 

• The metal-line statistics to be probed by COS rarely 
vary by a factor of more than two between the different feed- 
back models, despite enrichment patterns varying by much 
greater factors (e.g., 6.4 times more metals in the z — 
WHIM of cw vs. vzw). The relatively restricted density- 
temperature phase space traced by Ovi, Civ, Siiv, and 
Nevm are enriched similarly in all wind models relative to 
the total distribution of metals outside galaxies (§3.2,5.1.2). 

• The observed cosmic ion density Q summed from ab- 
sorbers typically underestimates the total Q of an ion. The 
stochastic effects of high overdensities being probed rarely in 
a volume-averaged measurement plus a large fraction of an 
ion's cosmic density being contained in saturated lines make 
accurate determinations of f^Clv ant ^ ^Siiv more difficult. 
We estimate recovery fractions of 56-69% for J1q vj and 30- 
39% for f^Ng vin * n S/N=30 COS observations. These frac- 
tions may be sensitive to the quality of the spectra and the 
fitting method used to estimate the column densities (§5.2). 

• Exploring alternative ionisation backgrounds given 
the observed range of quasar emission slopes beyond the 
Lyman limit, we find that the slope of the Ovi column 
density distribution becomes shallower using a harder 
background. The ionisation background strength varies 
by a factor of 16 at the O VI photo-ionisation edge (8.4 
Rydberg), but even for the weakest field Ovi remains 
primarily photo- ionised. Slopes of a = 1.57 and harder 
predict more photo-ionised Ovi arising from halo gas (§6.1). 

While metal lines observed by COS are more sensitive to 
the processes regulating galaxy evolution, they are also more 
uncertain in their modelling and interpretation compared to 
the Lya forest. Our parameter exploration finds good fits to 
existing O VI statistics, but appears to over-produce C iv 
and possibly under-produce Nevm arising in hot haloes. 
Like OD09, we favour a momentum-conserved wind model 
with turbulent broadening, although we note that metal-line 
cooling in the presence of the photo-ionising background re- 
duces O VI by a third and makes the average absorber hot- 
ter (30,000 vs. 15,000 K); however our Ovi results are still 
consistent with the current data. We further advocate this 
model for its ability to reproduc e the galactic stellar mas s 
function below M = 5x 10 10 M Q (jOppenheimer et alj201Ch , 
which are the same galaxies that are primarily responsible 
for IGM enrichment (OD09). 

We double down on our prediction of a primarily photo- 
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ionised origin for most UV absorbers (OD09), but realise 
the risky nature of our gamble. Our model variations show 
just how sensitive metal-line statistics are to the shape 
of the ionisation background at poorly constrained ener- 
gies, the stiff behaviour of metal-line cooling, and a tur- 
bulence model that is surely over-simplified. How metals 
mix in haloes and the IGM is perhaps the most under- 
appreciated uncertainty for simulating such observations. 
Our scheme of launching ISM gas with no mixing may 
be the largest difference with other works with smoother 
meta l licities leading to hotter temperatures JCen fc Chisaril 
120111 ; ISmith et all l201ll ; iTepper-Garcia et alj 120111 ); how- 
ever, there is observational evidence and theoretical motiva- 
tion that sma ll, highly-enriched cloudlets rema i n un- mixed 
in the IGM (|Simcoe et all 120061 ; ISchave et af1l2007T l. The 
other numerical uncertainty is the accuracy of the basic 
Gadget-2 method of implementing galactic winds by one- 
at-a-time particle ejection, also used in many other SPH 
codes. We have considered several variations of the scaling 
of outflow rate and ejection speed with galaxy mass, but we 
have not altered this basic scheme. 

The work presented here represents a broad brush ap- 
proach to simulating metal line absorbers along random 
sight lines sampling the volume-weighted Universe. With 
COS, new statistics can be generated including correlation 
functions and line of sight variance statistics, probing the 
stochasticity of various ion species and perhaps the phase of 
gas that they trace. Relating the distribution of metal ab- 
sorbers to the surrounding large-scale structure, particularly 
the low-mass galaxies responsible for enriching the IGM, will 
provide independent constraints on the nature of outflows. 
Studying the galaxy- absorber connection uniquely links the 
physics regulating the growth of galaxies to the dominant 
reservoir of primordial and enriched baryons in the IGM, 
which we will explore in the next paper in this series. 

We expect confronting the galaxy-absorber connection 
in particular with these and future simulations will provide 
unique tests that can both better distinguish wind models 
from each other, and push the limits of our scheme. Our 
claim that photo-ionisation can explain the majority of the 
high ion species explored here will face a significant test 
as we focus on the circumgalactic environment within sev- 
eral hundred kpc of galaxies. Given that interface physics 
and chemistry is not resolved in our simulation while it ap- 
pears fundamental for absorption in Galactic high-velocity 
clouds, suggests we are underestimating this phase of metal- 
absorbing gas. Conversely, this phase may represent a small 
minority of the universal quantity of diffuse metal-enriched 
absorbing gas, which is the view we advocate here. 
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